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were  computed  for  Lg  and  Pmax  for  earthquakes  and  explosions  in  the  EUS  and  WUS 
and  a  separation  of  0.6  magnitude  units  was  observed  thus  forming  a  regional 
discriminant.  The  Lg  is  larger  for  earthquakes  than  for  explosions.  This  con¬ 
clusion  is  somewhat  uncertain  in  the  WUS  because  of  the  large  scatter  in  the  ex¬ 
plosion  population.  However,  a  reasonable  explanation  for  this  large  scatter  is 
that  the  small  events  at  NTS  are  at  such  shallow  depths  in  dry  alluvium  that  the 
medium  is  weak,  resulting  in  a  low  corner  frequency.  This  decreases  the  ratio 
Pg/L  since  P  contains  comparatively  more  high  frequency  than  Lg.  The  scatter 
probably  would  not  be  a  problem  in  a  true  test  ban  situation  since  shots  will  be 
well  buried  to  avoid  surface  collapse. 

The  WUS  earthquake  P  versus  L  is  displaced  about  0.2  magnitude  units 
from  the  EUS  curve,  reflecting  the  lalge  P  phase  relative  to  L  in  the  WUS. 

The  GNOME  explosion  and  the  Hebgen  Lale  earthquake,  which  Ire  on  the  border 
between  the  EUS  and  WUS  defined  by  Der,  Mass£  and  Gurski  (1975),  show  differ¬ 
ences  in  amplitude  distance  relations  for  the  same  event  in  different  provinces. 
The  greater  WUS  attenuation  results  in  amplitudes  at  10°  of  about  0.4  magnitude 
units  below  that  in  EUS. 

Analysis  of  the  SALMON  and  18  February  1964  Alabama  earthquake  shows  that 
there  is  no  earthquake/explosion  discrimination  capability  using  maximum  trans¬ 
verse  to  maximum  radial  amplitude  ratios.  We  also  find  that  the  source  spectra 
of  the  two  events  are  identical  between  1  and  10  Hz,  that  the  L  spectrum  is 
different  from  the  P  spectrum  and  is  therefore  not  the  source  spectrum,  and  that 
the  L  spectrum  is  contaminated  by  scattered  coda  from  earlier  phases  so  that 
high-?requencies  observed  in  the  L  phase  may  not  be  predictable  by  any  deter¬ 
ministic  theory  of  L  .  Identical  lonclusions  with  respect  to  the  P  and  L 
spectra  were  obtainel  by  analysis  of  spectra  of  Soviet  explosions  as  observed  at 
NORSAR. 

In  the  EUS,  10  Hz  energy  is  observed  out  to  10°,  and  for  Soviet  shots  ob¬ 
served  at  NORSAR,  out  to  33°.  The  peak  in  the  S/N  varies  from  5  Hz  at  11°  to 
2  Hz  at  33°;  however,  for  small  events  which  have  not  yet  been  studied,  the  peak 
may  well  be  at  higher  frequencies  because  of  higher  corner  frequencies  and  be¬ 
cause  much  of  the  noise  on  LRSM  and  NORSAR  systems  at  10  Hz  may  be  system  or 
quantization  noise  which  can  be  reduced  by  more  carefully  designed  systems. 

Detection  thresholds  for  a  single  element  at  the  C3  subarray  of  NORSAR  were 
determined  to  be  m^  0.9,  2.4,  3.3,  and  3.9  at  5°,  10°,  15°  and  20°,  respectively. 
It  should  be  notedthat  the  lower  frequency  noise  levels  are  high  at  NORSAR, 
which  is  near  the  sea,  and  thus  that  stations  on  continents  might  have  lower 
thresholds  at  large  distances  where  the  lower  frequencies  are  the  most  useful  for 
detection. 
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ABSTRACT 

Amplitude-distance  curves  are  different  in  the  EUS  and  WUS:  P  decays 

-25-30  max 

as  r  '  and  r  *  in  EUS  and  WUS,  while  for  the  maximum  after  3.6  km/sec  on 

-2  -3 

the  vertical  component  (termed  L  )  the  decay  rates  are  r  and  r  .  The  EUS 

s 

results  are  in  general  agreement  with  the  literature  and  with  the  data  presented 
by  Nersesov  and  Rautian  (1964)  for  events  on  the  northern  margin  of  tectonic 
regions  in  the  Soviet  Union  suggesting  that  discrimination  results  in  the  EUS 
are  relevant  to  NSS  stations  within  the  Soviet  Union. 

Using  these  distance  amplitude  relations,  network  mean  amplitudes  at 

1000  km  were  computed  for  L  and  for  earthquakes  and  explosions  in  the 

EUS  and  WUS  and  a  separation  of  0.6  magnitude  units  was  observed  thus  forming 

a  regional  discriminant.  The  L  is  larger  for  earthquakes  than  for  explosions. 

8 

This  conclusion  is  somewhat  uncertain  in  the  WUS  because  of  the  large  scatter 

in  the  explosion  population.  However,  a  reasonable  explanation  for  this  large 

scatter  is  that  the  small  events  at  NTS  are  at  such  shallow  depths  in  dry 

alluvium  that  the  medium  is  weak,  resulting  in  a  low  corner  frequency.  This 

decreases  the  ratio  P  /L  since  P  contains  comparatively  more  high  frequency 

SB  £ 

than  L^.  The  scatter  probably  would  not  be  a  problem  in  a  true  test  ban  situ¬ 
ation  since  shots  will  be  well  buried  to  avoid  surface  collapse. 

The  WUS  earthquake  versus  is  displaced  about  0.2  magnitude  units 

from  the  EUS  curve,  reflecting  the  large  P^  phase  relative  to  in  the  WUS. 

The  GNOME  explosion  and  the  Hebgen  Lake  earthquake,  which  are  on  the 
border  between  the  EUS  and  WUS  defined  by  Der,  Mass£  and  Gurski  (1975),  show 
differences  in  amplitude  distance  relations  for  the  same  event  in  different 
provinces.  The  greater  WUS  attenuation  results  in  amplitudes  at  10°  of  about 
0.4  magnitude  units  below  that  in  EUS. 

Analysis  of  the  SALMON  and  18  February  1964  Alabama  earthquake  shows  that 
there  is  no  earthquake/explosion  discrimination  capability  using  maximum  trans¬ 
verse  to  maximum  radial  amplitude  ratios.  We  also  find  that  the  source 
spectra  of  the  two  events  are  identical  between  1  and  10  Hz,  that  the  L 

8 

spectrum  is  different  from  the  P  spectrum  and  is  therefore  not  the  source 
spectrum,  and  that  the  L  spectrum  is  contaminated  by  scattered  coda  from 
earlier  phases  so  that  high-frequencies  observed  in  the  L^  phase  may  not  be 
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predictable  by  any  deterministic  theory  of  L  .  Identical 

S 

respect  to  the  P  and  spectra  were  obtained  by  analysis 
Soviet  explosions  as  observed  at  NORSAR. 


conclusions  with 
of  spectra  of 


In  the  EUS,  10  Hz  energy  is  observed  out  to  10°,  and  for  Soviet  shots 
observed  at  NORSAR,  out  to  33°.  The  peak  in  the  S/N  varies  from  5  Hz  at 
11°  to  2  Hz  at  33°;  however,  for  small  events  which  have  not  yet  been  studied, 
the  peak  may  well  be  at  higher  frequencies  because  of  higher  corner  frequen¬ 
cies  and  because  much  of  the  noise  on  LRSM  and  NORSAR  systems  at  10  Hz  may 
be  system  or  quantization  noise  which  can  be  reduced  by  more  carefully 
designed  systems. 


Detection  thresholds  for  a  single  element  at  the  C3  subarray  of  NORSAR 
were  determined  to  be  m^  0.9,  2.4,  3.3,  and  3.9  at  5°,  10°,  15°  and  20°, 
respectively.  It  should  be  noted  that  the  lower  frequency  noise  levels  are 
high  at  NORSAR,  which  is  near  the  sea,  and  thus  that  stations  on  continents 
might  have  lower  thresholds  at  large  distances  where  the  lower  frequencies 
are  the  most  useful  for  detection. 
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INTRODUCTION 


In  1958  the  Geneva  Conference  of  Experts  suggested  that  a  worldwide 
network  of  180  seismic  stations  might  constitute  a  feasible  seismic 
monitoring  system  to  enforce  a  comprehensive  nuclear  test  ban  (CTB).  Such 
a  large  number  of  stations  naturally  included  several  inside  the  United 
States  and  the  USSR,  and  there  was,  therefore,  considerable  research  carried 
out  in  succeeding  years  on  distance-amplitude  relations  and  discrimination 
capabilities  at  "regional"  distances  of  less  than  say,  20°. 


In  time  it  became  less  clear  that  internal  stations  within  the  USSR  would 
be  allowed,  national  means  of  verification  were  emphasized,  and  furthermore, 
seismologists  realized  that  teleseismic  signals  were  not  so  dominated  by  effects 
of  complicated  crustal  structures  as  were  the  regional  signals.  Since  the 
teleseismic  signals  were  easier  to  understand,  they  could  perhaps  be  relied 
upon  more  for  discrimination.  Thus,  around  1963,  emphasis  shifted  from  studies 
at  regional  distances  to  those  at  teleseismic  distances. 

In  1977  it  began  again  to  seem  possible  that  the  USSR  would  allow  stations 
within  its  borders,  and  so  interest  has  been  revived  in  regional  discrimination. 
In  retrospect,  it  seems  clear  that  research  was  on  the  verge  of  major  discov¬ 
eries  on  the  subject  of  regional  discrimination  in  1963,  when  emphasis  was 
shifted  to  the  teleseismic  distances.  In  the  interim,  due  to  work  related  to 
earthquake  risk,  plate  tectonics,  and  to  some  degree  the  ARPA  discrimination 
program,  substantial  improvement  in  our  understanding  of  the  propagation  of 
crustal  phases  has  nonetheless  occurred. 


-This  paper  begins  with  a  review  of  the  literature.  Such  a  review  is  useful 
for  all  workers  in  the  field  as  it  aids  in  avoiding  repetition  of  work,  and 
shows  the  rather  advanced  state  of  existing  knowledge.  This  advanced  state  was 
not  recognized  by  the  author,  and  most  other  workers  in  the  field,  when  the 
subject  of  regional  propagation  and  discrimination  again  opened  up  in  1977. 


The  review  begins  with  work  on  propagation  of  P^,  Pg,  S^,  L^,  and  R^  and 

continues  with  the  history  of  discrimination  at  regional  distances  up  to  1978. 

New  work  performed  especially  for  this  report  is  discussed  next;  distance 

amplitude  relations  for  P^,  P^,  and  for  the  maximum  amplitude  observed 

before  and  after  S  are  presented  for  earthquakes  and  explosions  in  both  the 
n 

Eastern  United  States  (EUS)  and  the  Western  United  States  (WUS) .  Considerable 
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care  was  taken  in  these  analyses  to  avoid  bias  due  to  clipping,  to  discarding 
of  clipped  readings,  or  to  missing  the  maximum  excursion  during  examination 
of  a  film  record.  Use  was  made  of  the  low  gain  analog  tape  data  as  necessary. 
Bias  effects  are  maximized  at  close  distances  for  large  events  and  lead  to 
estimates  of  the  rate  of  decay  of  amplitude  with  distances  which  are  too  small. 
Care  was  also  taken  to  plot  noise  levels  on  distance  amplitude  plots  even  if  a 
signal  could  not  be  detected.  This  helps  avoid  the  additional  bias  toward 
a  small  rate  of  decay  which  can  occur  when  small  amplitudes  at  large  distances 
are  not  plotted  because  they  are  not  detected,  whereas  the  unusually  large 
amplitudes  are  detected  and  are  plotted. 

By  means  of  references  to  the  work  reviewed  in  the  literature  review, 
these  results  are  then  placed  in  context,  and  decay  rates  are  attached  to 
the  amplitude-distance  curves  so  that  they  may  be  used  in  magnitude  estima¬ 
tion.  Then  these  distance-amplitude  relations  are  used  to  define  event  magni¬ 
tudes  (log  amplitude  at  1000  km)  for  each  phase  and  the  discrimination  capa¬ 
bility  of  a  compressional  and  shear  magnitude  is  investigated  in  a  fashion 
analogous  to  the  conventional  M^m^. 

Finally,  several  topics  are  discussed  with  the  aid  of  spectra  of  explo¬ 
sions  and  earthquakes.  Among  these  are  detection  threshold,  optimum  filter  for 
detection  of  regional  events,  and  difference  in  source  spectra  for  earthquakes 
and  explosions. 


REVIEW  OK  LITERATURE 


Amplitude-Distance  and  Travel-Time 

Press  and  Ewing  (1952)  apparently  were  the  first  to  discuss  the  phases 
L  and  R  .  They  discovered  them  on  newly  installed  seismographs  at  Palisades, 

O  O 

New  Jersey,  and  associated  them  with  earthquakes  in  California.  The  instru¬ 
ment  response  was  nearly  flat  in  the  period  range  1  to  18  seconds,  and  they 
reported  dominant  periods  ranging  from  1/2  to  6  seconds  for  the  L  phase  which 
exhibited  primarily  transverse  motion.  The  signal  also  exhibited  reverse  dis¬ 
persion  with  the  shorter  periods  arriving  first.  The  arrival  of  these  short 
periods  was  quite  sharp,  with  a  group  velocity  3.51  ±  .07  km/sec. 

The  Rg  phase  arrived  at  3.05  ±  .07  kra/sec  with  periods  between  8  and  12 
seconds  at  maximum  amplitude.  It  exhibited  predominantly  radial  and  vertical 
motion. 


Another  experimental  investigation  was  reported  by  Blth  (1955).  He  used  a 

Wiechert  instrument  which  he  stated  was  not  very  sensitive  to  periods  of  less 

than  3.0  seconds.  He  determined  that  there  were  two  statistically  distinct 

L  phases,  arrivals  at  Uppsala  from  Eurasia;  he  referred  to  them  as  LI  and  L2. 

8  8  8 
The  mean  period  for  waves  was  LI: 5. 8  sec,  L2:6.8  sec,  R  :9.2  sec;  therefore, 

g  8  g 

these  observations  cannot  be  of  direct  interest  in  the  discrimination  problem, 

lying  as  close  as  they  do  to  the  peak  in  the  microseism  band.  The  velocity 

ranges  defined  by  B&th  for  the  various  phases  an:  Ll:  v  >  3.46  km/sec;  L2 

3.25  <  v  <  3.46;  R  :  v  <  3.25.  BSth  further  divided  his  population  of  L2 
8  8 
arrival  velocities  into  two  groups,  L2'  and  L2",  using  the  velocity  3.36  km/sec. 

8  8 

From  reading  his  paper  it  is  not  clear  to  the  present  authors  whether  these 
distinct  L^  arrivals  are  due  to  multipathing,  to  different  Airy  phases  due  to 
group  velocity  extrema  or,  in  some  cases  possibly  due  to  different  structures 
along  different  paths  to  Uppsala.  A  considerable  amount  of  plotted  velocity 
scatter,  as  contrasted  to  multiple  arrivals  observed  on  single  records  is  cer¬ 
tainly  due  to  inaccurate  epicenters. 

B&th's  median  L  velocity  was  3.45  km/sec  as  compared  to  3.51  given  by 
Press  and  Ewing.  Since  BSth  suggests  that  errors  of  0.15  km/sec  would  not  be 
unexpected  due  to  epicentral  errors,  the  two  results  seem  to  be  within  experi¬ 
mental  error  of  each  other.  BSth's  median  R^  velocity  is  3.07,  in  good  agreement 
with  the  3.05  value  of  Press  and  Ewing. 
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Press  and  Gutenberg  (1956)  searched  records  from  the  Kern  County  earthquake 
of  July  21,  1952  for  the  channel  P  wave  phase  analogous  to  the  channel  S  waves 
which  they  hypothesized  to  make  up  L  .  They  were  successful,  finding  a  phase 
with  a  mean  velocity  of  6.09  km/sec  out  to  a  distance  of  37.5°.  They  suggested 
the  name  irg.  In  general  they  characterized  it  as  having  small  amplitudes.  Only 
two  of  their  observations  were  at  distances  of  less  than  12°. 

I.  Lehman  (1957)  studied  the  properties  of  L^  as  read  on  North  American 
records.  Hers  appears  to  be  the  first  discussion  of  the  signal  character  as 
seen  on  Benioff  short-period  instruments.  She  observed  that  signals  recorded  at 
Palisades  from  an  Oklahoma  event  had  periods  of  1/2  to  1  second.  Observations 
in  Pasadena  of  the  same  event  had  longer  periods  and  apparently  smaller  amplitudes, 
although  calibrations  were  not  completely  reliable.  She  speculated  that  this  might 
be  related  to  the  larger  felt  areas  of  earthquakes  in  the  East. 

Oliver  and  Ewing  (1957,  1958a, b)  discussed  the  foregoing  observations  of  L^, 
and  some  of  their  own,  in  terms  of  fundamental  Love  wave  dispersion  and  of  funda¬ 
mental  and  higher  mode  Rayleigh  dispersion.  They  identified  the  LI  and  L2  of  BSth 
with  successive  extrema  at  5  and  6  seconds  period  of  the  first  higher  Rayleigh 
mode;  and  pointed  out  that  it  would  be  only  a  coincidence  of  the  structure  that  the 
velocities  at  these  extrema  were  near  the  crustal  shear  velocity  to  which  this 
higher  Rayleigh  mode  tends  in  the  high-frequency  limit  and  which  is  the  Lg  velocity 
for  high  frequencies.  R^  for  periods  greater  than  2  seconds  they  identified  with 
a  broad  plateau  in  the  fundamental  dispersion  diagram;  they  showed  how  sediments 
could  greatly  reduce  the  fundamental  Rayleigh  mode  velocities  at  high  frequencies, 
resulting  in  a  long,  dispersed  coda. 

Romney  (1959)  made  use  of  VELA  Uniform  data  to  report  some  of  the  first 
reliable  amplitude  distance  data  at  regional  distances.  The  data  were  taken 
along  a  profile  extending  from  the  Nevada  Test  Site  (NTS)  Southeast  into  Texas 
and  thence  Northeast  to  Maine.  In  Figure  1  we  see  the  data  (for  the  phase  Pn) 
and  note  that  out  to  800  km  a  line  drawn  as  the  inverse  cube  of  the  distance 
provides  a  good  fit. 

Romney  also  discussed  the  L  waves,  noting  that  they  propagated  with  a 
velocity  of  about  3.5  km/sec.  In  Figure  1  we  see  the  amplitudes  of  these  waves 
on  all  three  components.  Again,  an  amplitude  decay  as  the  inverse  cube  is 
suggested  by  Romney  as  appropriate  out  to  the  limit  of  observation  (-2.5  fits 
better);  and  all  three  amplitude  components  seem  equal  to  one  another. 
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-Amplitude  of  P  waves  from  a  19-kt  un¬ 
derground  explosion  as  a  function  of  distance.  The 
dashed  curve  is  drawn  proportional  to  the  inverse 
i-ulie  of  the  distance. 


Amplitudes  of  IJWcm/stc  abaar  waves  from  a  19-kt  underground  explosion.  Curves 
proportional  to  the  inverse  cube  of  the  distance. 


Figure  1.  Distance-amplitude  relations  from  Romney  (1959)  for  P  and  L 
waves .  ® 


Shurbet  (1960)  analyzed  the  irg  (Shurbet  preferred  to  call  it  P)  phases  of 
Press  and  Gutenberg  (1956)  further  and  concluded  that  the  reason  they  did  not 
propagate  strongly  to  the  east  was  the  thinning  of  the  crust.  He  relied  on 
work  by  Ewing  and  Press  (1959)  stating  that  the  M  discontinuity  is  at  47  km 

beneath  the  Rockies  and  shallows  to  35  km  beneath  the  interior  plains.  If  the 

attenuation  of  a  phase  is  strongly  controlled  by  structure,  it  clearly  becomes 
important  to  be  very  careful  about  distance-amplitude  relations. 

Romney  et  al  (1962)  published  substantial  travel-time  and  amplitude  data 
for  the  event  GNOME.  In  Figures  2a-e  their  data  shows  how  the  P  travel  times 
could  be  contoured  at  regional  distances,  and  how  the  trend  of  the  P  (P)  travel- 
times  paralleled  the  P  times  and  seemed  to  have  only  slightly  greater  variance 
about  a  smooth  line  drawn  through  the  data  than  does  P^.  With  L^,  as  pointed 
out  by  Romney  et  al,  the  situation  is  different.  There  seems  to  be  much  greater 
scatter  and  no  apparent  correlation  with  P^  and  P  travel  times.  Also  in  Figures 

2d  and  e  we  reproduce  the  amplitude-distance  data  presented  by  Romney  et  al. 

Brune  and  Dorman  (1963)  used  signals  from  earthquakes  in  Northern  Canada 
as  recorded  in  Southeast  Canada  and  at  Palisades  to  extend  the  work  of  Oliver 
and  Ewing  discussed  above,  and  also  to  investigate  the  properties  of  Pn  and  Sn> 
Figure  3  gives  some  of  the  relevant  data. 

Ryall  and  Stuart  (1963)  measured  arrival  times  and  amplitudes  from  eight 
NTS  explosions  along  a  profile  to  Ordway,  Colorado.  A  short  spread  of  instru¬ 
ments  was  set  out  at  28  locations  along  the  profile  so  that  phase  velocity  could 
be  measured.  Ryall  and  Stuart’s  distance-amplitude  data  for  P  and  P  are  repro- 

n  _3  _4 

duced  in  Figure  4.  While  Romney  et  al  (1962)  apparently  favored  r  and  r  for 
these  two  phases,  best  fits  to  Ryall  and  Stuart's  data  (on  a  different  profile, 
of  course)  yield  r  -3.28  and  r  -3.52. 


Press  (1964)  analyzed  LRSM  recordings  of  P^  and  as  detailed  in  Figure  5. 

In  deriving  his  Q  values  he  assumed  for  both  P  and  L  that  the  amplitude  decayed 

in  the  time  domain  as  r  exp(-irr/UQT)  with  U  =  3.5.  (This  is  clearly  incorrect 

for  P  for  which  U  should  be  about  6.0.) 
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I  M  l  ^  It  si  \  I  I  S 


Travel-time  and  distance  amplitude  relations  for  the  event 
GNOME  from  Romney  et  al  (1962). 


DISTANCE  (KM) 


Comparison  o(  Pn  and  P  travel  time  residuals  relative  Vo  constant  velocities  (or 
profile  east 


Figure  2b.  Travel-time  and  distance  amplitude  relations  for  the  event 
GNOME  from  Romney  et  al  (1962) . 
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Kirkland  Lake  seismograms  of  t lie  Baffin  Bay  earthquake  of  May  2,  1957. 
Tlie  epicenter  is  nearly  due  nortli  of  Kirkland  Lake. 


Travel  time  data  for  the  Canadian  shield. 

Figure  3.  Signal  characteristics  and  travel  times  for  Pr  and  Sn  from 
Brune  and  Dorman  (1963) . 
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I )»«  rtrMilii*  ol  :t  Kt  tidy  of  .si  iami'-  ttavrs  genera  UmI  bv  eight  noi'lrar  explosions  and 
recorded  at  31  locations  U'tween  the  Nevada  teat  site  (NTS)  und  Ordwa\ ,  Colorado,  are 
discussed  The  line  of  recording  stations  crosses  the  eastern  part  of  the  Basin  and  Range 
province,  the  Colorado  plateaus,  and  the  southern  Rocky  Mountains,  and  it  extends  into  the 
Great  Plains  In  the  eastern  Basin  and  Range  province  and  the  western  part  of  the  Colorado 
plateaus  (0  <  A  <  385  km),  the  time-distance  curves  for  P,  and  Pm  can  be  expressed,  respec¬ 
tively,  as  Ti  ~  08  -f-  A/6  0  and  T*  =  58  +  A/7  6  A  third  phase,  tentatively  identified  as 
P*.  is  lepreaented  by  the  equation  T»  =  38  4-  A/65.  Using  the  crustal  structure  and  Pm  veloc¬ 
ity  (79  km/sec)  found  for  the  NTS  region  by  other  authors,  we  find  that  the  above  relations 
indicate  that  the  thickness  of  the  crust  increases  from  about  25  km  at  NTS  to  about  42  km 
in  the  western  part  of  the  Colorado  plateaus.  East  of  this  boundary  the  velocity  of  P  in  the 
upper  mantle  increase  to  SO  km /sec ;  depth  to  the  M  discontinuity  in  the  Colorado  plateaus 
is  approximately  constant  over  the  range  435  <  A  <  645  km.  Beyond  850  km,  first  arrivals 
in  the  southern  Rocky  Mountains  ami  Great  Plains  provinces  indicate  an  apparent  velocitv 
of  about  8  4  km/sec  Amplitudes  of  Pm  attenuate  according  to  the  equation  A  =  A«A-,/e 
(A  —  d)~M  e~*  over  the  distance  range  150  C  A  <  850  km  (d  =  GO  km)  This  relation 
yield*  a  value  of  <?,  for  Pm,  of  about  520  The  amplitudes  of  P ,  attenuate  extremely  rapidly, 
and  beyond  about  130  km  this  phase  cannot  be  identified  with  certainty.  An  extension  of  the 
P,  travel-time  branch  at  long  distances  could  be  associated  with  svaves  reflected  beyond  the 
critical  angle,  from  the  bane  of  the  crust  This  phase,  railed  P  nfter  Mohornvific,  appears  io 
attenuate  as  A  —  A%  c~*mr%A A'1'*  The  value  of  Q  indented  bv  this  equation  is  about  200 


Nuclear  Test  Data 


Shot 

Name 

Date 

(1962) 

Location  Origin  Time,  UT 

Latitude  Longitude 

N  W  h  rn  s 

Yield, 
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M  edium 

I 
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00 
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4  5 
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11 
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30 
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I  8 
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Cimarron 

Feb.  23 

37*07.7'  1 16*02.9'  18 

00 
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IV 
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■I 
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Estimated  Frequencies  (/,  cpsj  and  Measured  Amplitudes  (A,  m*i) 

Amplitudes  for  shots  I  to  V  were  scaled  to  the  Hardhat  explosion  using  shot  factors;  amplitudes 
for  shot  VI  were  shifted  by  an  arbitrary  amount  to  approximately  scale  to  Hardhat 
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3  0 
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Figure  4a.  Signal  characteristics  and  amplitude  distance  relations  from 
Ryall  and  Stuart  (1963). 
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Figure  4b.  Signal  characteristics  and  amplitude  distance  relations  from 
Ryall  and  Stuart  (1963). 
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Figure  5.  Q  and  the  implied  amplitude  distance  relationships  for  P 

L  for  various  profiles  from  Press  (1964). 
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Press'  SE  profile  has  four  events,  and  Q  varies  from  370  to  510  for  L 

8 

and  from  134  to  340  for  P  .  Apparently  there  is  a  fair  amount  of  scatter 
inherent  in  determinations  of  decay  rates.  The  lowest  Q  of  134  corresponds 

to  a  period  of  0.82  yielding  a  QT  product  of  109.  This  may  be  converted  into 

-4.6  — 

a  decay  rate  r  '  which  would  fit  fairly  well  the  observations  for  P  obtained 

from  GNOME  by  Romney  et  al  as  seen  in  our  Figure  2. 

Nersesov  and  Rautian  (1964)  presented  comprehensive  results  from  a  string 
of  stations  along  the  southern  border  of  the  USSR  substantially  within  the 
seismically  active  region  of  the  nation.  The  profile  extended  from  the  Pamirs 
on  the  Southwest  to  Lake  Baikal  on  the  Northeast. 

An  event  on  January  31,  1962  (at  00:05:57,  38.5°N,  70.3°E  at  60  km  depth, 

no  reported  magnitude  according  to  NEIS)  served  as  a  reference  event  for  this 

study.  In  Figures  6a-d  we  see  amplitudes  of  S(L  ),  P(P  ),  P  and  P  for  this 

8  8  n 

and  other  events  as  a  function  of  distance.  Analysis  of  the  L  data  gives  a 

-2.1  -2.3  8 

decay  rate  of  r  ‘  out  to  3000  km,  and  r  *  out  to  3500  km.  This  is  in 

good  agreement  with  average  results  derived  from  Press  (1964) ,  but  is  a  smaller 

_3 

rate  of  decay  than  the  r  suggested  by  Romney  (1959) . 

In  comparison  wit  nney  et  al  (1962)  the  decay  rate  is  much  less  than 
-4 

that  of  r  found  on  the  Western  profile  out  of  GNOME,  but  is  in  agreement 
with  the  Eastern  profile. 

Calculations  based  on  Nersesov' s  composite  amplitude  graphs  give  as  decay 
rates:  L^,  2.28  to  3000  km;  P  ,  2.43  to  1500  km.  Distance  amplitude  averages 
for  Pn  and  P  seem  to  be  too  erratic  for  meaningful  averaging. 

Nersesov  noted  that  the  maximum  compressional  amplitude  on  the  record 
seemed  to  be  a  more  stable  measure  of  magnitude  than  any  particular  phase, 
amplitude  is  plotted  in  Figure  6d  as  two  lines,  one  for  earthquakes  from 
Baikal  region  and  one  from  elsewhere.  He  also  has  plotted  a  maximum  shear 
wave  amplitude.  To  a  distance  of  2000  km  the  slope  for  the  shear  wave  is  2.24, 
and  the  slope  for  the  compressional  wave  maximum  is  2.3;  with  the  compressional 
wave  having  perhaps  0.6  logarithmic  units  less  amplitude  than  the  shear  maximum. 
(We  shall  see  that  for  explosions  the  shear  amplitude  would  presumably  be  more 
comparable  to  the  compressional  amplitude.) 

Haskell  (1964)  attempted  to  explain  theoretically  the  P  (P)  amplitude 
decay  rates  found  by  Press  (1964).  He  used  elementary  ray-theoretical  methods 
to  build  up  a  picture  of  P^  as  a  superposition  of  leaking  compressional  modes. 
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Composite  average  travel  time  curve.  1-West;  2-Siberia; 
numbers  on  the  curve  are  apparent  velocity  (km/sec) 

Figure  6a.  Travel-time  and  amplitude-distance  figures  from  Nersesov  and 
Rautian  (1964). 


Composite  amplitude  curves  for  earthquakes  of  the  Altai  Sayan 
recorded  to  the  southwest  (a)  .  and  to  the  northeast  (b). 


Composite  amplitude  curves  for  earthquakes  of  Pribaikal. 

Figure  6c.  Travel-time  and  amplitude-distance  figures  from  Nersesov  and 
Rautian  (1964). 


He  found  that  the  existence  of  P  could  be  associated  with  the  occurrence  of  a 
common  maximum  in  the  group  velocity  curves  for  high-order  multiple-P-ref lec¬ 
tion  modes.  It  is  notable  that  Haskell  attempted  to  explain  the  Q  quoted  by 
Press  of  Q  =  260;  but  since  Press  used  too  low  a  group  velocity  for  P  ,  the 
correct  value  for  Q  requiring  explanation  should  be  Q  =  260  x  (3. 5/6.0)  =  152. 
Haskell's  theory  was  ray-theoretic;  his  calculated  absorption  depended  exponen¬ 
tially  upon  the  number  of  free-suitace  reflections  encountered  by  each  ray. 

Thus,  the  absorption  was  independent  of  frequency,  and  therefore,  the  estimated 
Q  increased  linearly  with  frequency.  Haskell  estimated  a  Q  of  71  at  1  Hz; 
however,  the  observed  dominant  frequency  in  Press'  data  was  1.4  Hz,  and  values 
of  2  Hz  are  typically  seen  to  be  dominant  in  spectra.  Thus  one  might  as  well 
have  asserted  that  the  theoretical  Q  was  140,  in  agreement  with  the  Q  (as 
corrected  in  this  study)  measured  by  Press.  Thus,  Haskell  was  not  necessarily 
correct  in  asserting  that  a  low-velocity  surface  layer  is  necessary  to  bring 
Pg  theory  and  observation  into  agreement. 

The  situation  is  complicated  even  further  by  the  extreme  dependence  of 
absorption  on  group  velocity.  At  a  distance  of  1000  km  the  group  velocity  for 
which  Haskell  calculated  his  Q  value  arrives  2.8  seconds  after  the  limiting 
velocity  of  6  km/sec,  while  an  arrival  with  twice  the  Q  (due  to  fewer  reflec¬ 
tions)  arrives  after  1.4  seconds.  Overall  it  does  not  seem  possible  to  compare 
theory  and  observation  using  Haskell's  theory  except  in  the  sense  that  it 
predicts  a  sharp  arrival  near  the  limiting  velocity  of  P  in  the  upper  crustal 
layer.  One  must  then  attribute  the  long  P^  coda  often  observed  to  multipathing, 
or  to  phases  associated  with  multiple  shallower  layers  or  to  other  modes  such 
as  multiple  P  with  one  or  two  S  legs. 

Evernden  (1967)  published  travel-time  and  distance-amplitude  results  for 

explosions  in  the  WUS  and  EUS.  In  the  EUS,  he  plotted  P^  log(A/T)  versus  log 

distance  and  obtained  the  slopes  shown  in  Figures  7a-d.  The  average  slope  was 

-2.0 

r  ’  ,  and  his  magnitude  relation  was,  therefore,  as  given  in  Figure  7a  for 


Evernden  asserts  that  investigation  of  P  distance  amplitude  plots  in 

8  -3 

many  LRSM  shot  reports  indicates  that  P^  decays  in  the  WUS  as  r  .  Plotting 
the  interpolated  P^  amplitude  at  500  km  versus  m^  determined  from  regional 
observations  using  regional  distance-amplitude  formulas  he  produces  Figure  7c. 
Regression  on  yields  a  slope  of  1.14;  however,  considering  that  both  axes 
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m,w  «*  mt.i  *  -3.27  -f  log  A/T  +  2  log  A. 

Figure  7a.  Travel-time  and  amplitude-distance  figures  from  Evernden  (1967). 
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Figure  7b.  Travel-time  and  amplitude-distance  figures  from  Evernden  (1967) 
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Pq  Al 


mg.o  =1.28  +  1.14  log  (A/T) 500  (9) 

mg.o  =  -8.00  +  1.14  log  A/T  +  3.42  log  A  (10) 

mg.Q  *  1.6  +  log  (A/T)  +  3  log  (A/500)  (Interpretation  of  this  paper) 


Earthquakes  (A/T)i—  mb' 


Colons  (5  Feb  1962) 

40 

3.1 

3.11 

Hebgen  (25  Feb  1962) 

50 

3.2 

3.22 

Cache  Creek  (30  Aug  1962) 

2500* 

4.3 

5.16 

Cache  Creek  A/S  (5  Sep  1962) 

500* 

4.1 

4.36 

Bridgeport  (5  Apr  1962) 

42 

3.03 

3.29 

Fallon  (20  Jul  1962) 

240 

4.4 

4.16 

•  Based  on  W0S  data. 


Figure  7c.  Travel-time  and  amplitude-distance  figures  from  Evernden  (1967) 
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have  error  we  cannot  reject  a  slope  of  unity.  Thus  in  place  of  Evernden's 
formula  (10)  as  seen  in  Figure  7c  we  would  write  the  expression  immediately 
below  it. 

In  the  WUS  Evernden  concluded  that  the  Pfl  phase  was  erratic,  giving 
different  amplitude-distance  curves  on  different  profiles,  as  Figure  7d 
attests.  Evernden  also  asserted  that  the  partition  of  amplitudes  between 
refractors  varied  according  to  source  parameters. 

In  1969  Molnar  and  Oliver  (1969)  outlined  the  regions  of  favorable  propa¬ 
gation  for  the  phase  S^.  In  Figure  8  from  Molnar  and  Oliver  propagation  across 
the  Basin  and  Range  Province  in  the  WUS  is  seen  to  be  non-existent.  There  is  an 
indication  that  propagation  is  good  further  North  along  the  West  Coast  of  North 
America,  and  then  becomes  poor  again  along  the  concave  side  of  the  Aleutian  Arc 
and  into  the  Fairbanks  region  of  Alaska.  Obviously,  however,  more  data  is  needed 
to  fill  out  this  picture. 

I  Propagation  across  the  Baltic  Shield  seems  excellent,  but  through  the 

mountain  ranges  of  Southern  Europe  and  Asia  it  is  poor  and,  in  particular, 

t 

)  Molnar  and  Oliver  state,  "Moreover,  paths  crossing  mountainous  regions  of  Iran, 

I 

I  Pakistan,  and  Turkey  do  not  transmit  S  efficiently.  There  are  some  examples 

I  n 

;  of  S  at  MSH  for  short  paths  from  the  Hindu  Kush,  but  in  general,  predominant 

n 

frequencies  at  MSH  are  lower  than  at  NDI  (on  the  Indian  Shield)  for  these 
earthquakes. " 

This  quote  would  suggest  that  the  SRO  station  MAIO  presently  co- located 
near  MSH  would  be  somewhat  useful  although  not  ideal  for  studies  of  earthquakes 
within  the  USSR. 

In  general,  Molnar  and  Oliver  discovered  that  S^  is  not  well  propagated 
through  regions  which,  on  other  grounds,  are  suspected  of  high  attenuation 
due  to  melting. 

Figure  9  shows  results  from  Bollinger  (1970)  giving  travel  times  for  Pn> 

P  ,  and  L  out  to  1000  km  in  the  EUS.  In  general,  the  P  amplitude  data  are 
too  sparse  to  permit  determination  of  an  amplitude-distance  relation. 

I  In  Figure  10  we  see  Gumper  and  Pomeroy's  (1970)  travel-time  curve  for 

Africa. 
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Earthquake  of  May  31,  1966 
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Earthquake  of  March  8,  1966 
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*  Aeadmp  obtained  from  U9CGS 


Amplitude /period  versus  eg 
earthquakes.  Solid  circles— May  81,  lUv, 

1068;  Cross— April  36, 1086;  X— October  28, 1088. 


picentral  distance  for  the  P§  phase  of  central  Appalachian 
088;  Haif-eolid  circle— Nov.  34,  1964;  Circle -dot— March  8, 


Figure  9.  Travel-time  and  amplitude-distance  relations  from  Bollinger  (1970) 


Tr-avel-urne  dala  for  the  African  continent.  The  travel-time  curves  T  =  a/3  70  and 
T=  13  -f  A  ,'4  717  for  the  Canadian  shield  are  shown  for  comparison 


Figure  10.  Travel-time  data  for  Africa  from  Gumper  and  Pomeroy  (1970). 
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Baker  (1970)  (see  Figure  lla-c)  has  published  the  most  comprehensive 

data  set  for  ;  his  data  were  taken  from  the  LRSM  shot  reports  of  Nevada 

Test  Site  nuclear  explosions.  He  concluded  that  the  magnitudes  determined 

from  L  showed  substantially  less  scatter  than  did  magnitudes  determined  from 

P  •  Baker  fitted  distance-amplitude  relations  to  his  event-normalized  data 

over  several  distance  ranges  as  seen  in  Figure  11a.  Perhaps,  the  most  trouble- 

free  distance  range  is  that  of  200  to  2000  km.  At  smaller  distances  the  number 

of  stations  is  small;  and  possibly,  some  of  the  data  could  be  clipped  or 

amplitude  measurements  might  be  low  because  the  maximum  excursion  was  too  faint 

on  the  film.  At  larger  distances  amplitudes  are  often  enhanced  by  noise,  only 

larger  amplitudes  are  sometimes  observed,  and  we  move  from  the  WHS  to  the  EUS 

where  there  are  different  amplitude-distance  relations.  In  the  200  to  2000 

km  distance  range  the  fitted  slope  is  -2.3,  in  good  agreement  with  results 

derived  from  Press  (1964)  and  Nersesov  and  Rautian  (1964).  Baker  also  presents 

magnitudes  based  on  L  ,  Figure  lib.  He  also  calculated  L  station  corrections 

8  S 

(Figure  11c)  and  concluded  that  there  is  no  apparent  geographical  pattern  to  them. 

After  Baker's  study  there  are  few  publications  in  this  field  until  Nuttli's 

(1973)  study  of  earthquakes  with  epicenters  near  St.  Louis.  Nuttli  chose  to 

analyze  the  L  phase  as  an  Airy  phase,  and  therefore,  fit  the  data  with  an 
8  -1/3  -1/2 

equation  of  the  form  A  (sin  A)  exp(-yA),  where  A  is  distance  in  degrees. 

He  concluded  that  the  data  could  not  be  fit  adequately  by  a  straight  line  in 
a  log(amplitude) : log(distance)  plot.  However,  since  such  expressions  are 
common  in  the  literature,  he  did  fit  an  expression  of  this  form  as  seen  in 
Figure  12.  Nuttli  was  the  first  author  to  attach  an  absolute  magnitude  to  his 
distance-amplitude  scale.  He  accomplished  this  by  connecting  the  magnitude 
to  that  of  P^  derived  using  the  EUS  formulas  of  Evernden  (1967). 

Nuttli's  slope  of  -1.66  for  was  not  derived,  but  was,  presumably, 
selected  by  analogy  with  the  teleseismic  formula  for  LR.  Examination  of  Nuttli's 
Figure  12,  shows  that  while  one  could  not  reject  the  slope  of  -1.66,  the  data 
do  admit  an  interpretation  of  at  least  -2.0  for  the  distance  range  3°  to  30°. 

Even  at  closer  distances,  the  only  deviation  from  a  slope  of  -2.0  plotted  is 
that  for  the  October  21,  1965  event  at  distances  near  1.5°.  One  would  be 
tempted  to  suspect  that  possibly  at  such  close  distances,  the  maximum  excur¬ 
sion  might  be  missed,  or  that  data  from  clipped  stations  existed  that  could 
not  be  included  in  the  plots. 
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Slope  Values  of  Log  A/T  vehsus  Log  r  fob  Several  Distance  Ranges  of  t 


Approximate  lUnft  o(  rfkra)  Slope  of  Best  Fit  Line  for  All  Date  Used*  Root  Me»n  Square  Error 


50-250 

-1.548 

0.26788 

250-500 

-3.459 

0.23219 

200-1100 

-2.710 

0.26978 

200-2000 

-2.300 

0.29621 

300-4000 

-2.162 

0.30601 

500-1500 

-2.444 

0.23944 

1100-2500 

-1.537 

0.27357 

1500-2500 

-2.628 

0.23501 

2000-4000 

-1.551 

0.31006 

2500-4000 

-1.707 

0.30473 

•  One  slope  and  several  intercepts— one  for  each  event— were  obtained  for  each  distance 

range. 


A  ■  f  Uml 

Plot  of  reported  log  A/T  values  versus  epicentral  distance. 


Q  =  -0.14  +  (6.90/10>)r  -  (6.56/10V 
+  (3.37/10V  -  (9.38/10“)^ 

+  (l^/lO1*)^  -  (7.74/10 !,)r*. 


M  =  log  (A/T)  +  Q  +  s 

Figure  11a.  Amplitude-distance  data,  event  magnitudes  and  station  corrections 
for  L  from  NTS  from  Baker  (1970). 
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Figure  lib.  Amplitude-distance 

data,  event 

magnitudes 

and  station  corrections 

for  L  from  NTS  from  Baker  (1970). 
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Average  Magnitudes  (Normalized  M  -  5.0)  and  Residuals  for  Stations  Recording 
Five  or  More  NTS  Events 
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14 
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0.11 

L C  NM 

Lae  Cnicee,  N  M. 

4.80 

84 
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0.117 
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0.161 
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Figure  11c.  Amplitude-distance  data,  event  magnitudes  and  station  corrections 

for  L  from  NTS  from  Baker  (1970). 
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Observed  A/T  values  of  the  Z  component  of  l-*ec-penod  L,  wave*  for  the  earth- 
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solid  line  is  a  theoretical  attenuation  curve  for  an  Airy  phase  with  y  =  007  deg  '.  The  dashed 
line  is  a  straight  line  approximation  to  the  theoretical  curve  over  the  distance  interval 
3*  <  A  <  30* ;  over  these  distances  it  does  not  differ  from  the  theoretical  curve  by  more 
than  about  0.1  magnitude  unit.  The  arrows  indicate  0.1  magnitude  unit. 


Body  Wave  and  Surface  Nave  Magnitude*  of  Earthquakes 


20-Sec  Talesaisaic  3-  to  12-Sec  Teleseisaic 

Data  Surface  Saves  Rayleigh  Naves  P  Naves-  l-$ec  1  Naves 

9 


Nov.  9.  1968  5.20  2  0.28  (2) 

Oct.  21.  1965  4.13  i  0.32  (9) 

All p.  M,  1965 
*!«r.  3.  I 


S .19  *  0.19  (12) 
4.09  t  0.20  (31) 
2. 54  (3) 

4 .08  ♦  0.20  (13) 


5. SO  t  0.40  (30) 
4.85  t  0.23  (18) 
3.81  ♦  0.29  (9) 
4.76  -  0.33  (6) 


5.3S  1  0.23  (8) 

S .04  •  0.12  (14) 
3.83  ;  0.39  (12) 
4.67  ♦  0.10  HO) 


NvwbeTs  in  parentheses  refer  to  mafcer  of  observations. 

v*lu#*  obtained  froa  Pn  data  in  the  eastern  United  States  by  using  Euemden's 

(1967)  fonaule. 


Dete 


Coaperison  of  Body  Neve  Magnitudes  Obtained  froa 
Lg>  pn >  and  Teleseisaic  P  Naves  for  Earthquakes 
in  the  Central  United  States 


Origin  Tiae, 
UT 


Location 


% _ 

Pn  Teleseisaic  P 


Sep. 

IS,  1972 

OSh 

22a 

15. S* 

41.6*N,  89  .3*N 

4 .59 

(7) 

4.39 

(7) 

June 

4.  196? 

16h 

14a 

12.6s 

33.S*N ,  90  .9*N 

4.4) 

(!) 

4.49 

P) 

Jan. 

1,  1969 

23h 

3Sa 

36.2s 

34  ,8*N ,  92.6*N 

4. SO 

(9) 

4.30 

(9) 

July 

21,  1967 

09h 

14a 

48s 

37.S*N,  90 .6*N 

4.3? 

(10) 

4.21 

(12! 

Oct. 

1,  1971 

18h 

49a 

38.7s 

3S.8*N,  90  .4*N 

4.11 

(*> 

4.14 

(6) 

Feb 

12,  1971 

12h 

44a 

2  7.2s 

38.5*N.  87 .9"N 

3.33 

(«) 

3.24 

(S) 

June 

9.  1972 

19h 

ISa 

19.  is 

37. 7*N ,  90. 4* N 

3.27 

(5) 

2.96 

O) 

June 

19,  197? 

OSh 

46a 

14.7* 

J7.0*N,  89. 1*N 

3.2? 

fS) 

2.96 

0) 

Nuaber*  in  parentheses  refer  to  the  nuaber  of  observations. 


Figure  12. 


Amplitude-distance  and  event  data  for  L 


from  Nuttll  (1963). 
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Bollinger  (1973)  published  data  points  for  L^  versus  distance  for  three 
earthquakes  and  concluded  that  he  could  not  reject  Nuttli's  slope  of  -1.66. 

However,  we  see  that  in  Figure  13  the  right-hand  graph  is  really  two  earth¬ 
quakes  plotted  as  one.  If  the  1969  earthquake  points  are  dropped  by  0.4 
magnitude  units  as  is  appropriate  according  to  Bollinger's  Table  VI  m^  values 
of  3.93  and  3.59  then  the  least  squares  slope  becomes  -2.03  instead  of  -1.63. 

In  addition,  one  could  not  reject  a  slope  of  -2.0  through  the  1969  W.  Virginia 
earthquake . 

Street  (1976)  determined  slopes  greater  than  1.66  for  four  earthquakes  in 
the  St.  Lawrence  area.  While  some  of  the  lowest  amplitude  points  are  from  WUS 
stations,  even  in  the  East  the  slope  is  greater  than  1.66  -  perhaps  approaching 
2.5. 

In  1977  Jones,  Long  and  McKee  published  a  study  of  attenuation  in  the 
Southeastern  United  States.  In  Figure  14  we  reproduce  their  list  of  South¬ 
eastern  United  States  earthquakes,  together  with  reproductions  of  the  absolute 
distance-amplitude  curves  of  Richter  and  Nuttli.  It  seems  to  the  author  that 
one  could  not  reject  a  slope  of  -2.0  between  4°  and  9°.  However,  at  4°  there 
is  a  striking  discontinuity  in  the  data  which  the  authors  do  not  discuss  in 
detail;  between  4  and  5  degrees  the  amplitude  appears  to  increase  by  0.5 
magnitude  units.  This  seems  to  be  unusual  behavior  for  a  guided  or  surface  wave. 

The  data  consists  of  201  observations  from  72  events;  less  than  3  observations 

per  event.  It  seems  to  the  present  authors  that  this  few  observations  per  event 

could  well  lead  to  several  biases  of  a  statistical  nature.  Those  of  clipping, 

signal  enhancement  by  noise,  and  missing  small  observations  at  large  distances 

have  already  been  discussed.  The  nature  of  errors,  if  any,  which  might  result 

when  only  two  observations  are  available  (one  of  which  must  be  "used  up"  in 

computation  of  the  event  magnitude)  are  unknown  to  the  present  authors,  but  the  s 

i 

procedure  seems  intuitively  somewhat  unstable.  • 

t 

Ruzaikin  et  al  (1977)  have  published  a  valuable  survey  of  the  propagation  • 

characteristics  of  L^  in  Asia,  using  a  technique  very  similar  to  that  used  by  j  i 

Molnar  and  Oliver  (1969)  for  S^.  As  we  shall  see  later  in  this  report,  the 

amplitude  ratio  of  the  maximum  motion  after  Sr  to  the  maximum  amplitude  before 

S  is  a  powerful  regional  discriminant,  presumably  because  it  is  a  measure  of 
n 

the  ratio  of  shear  to  compressional  energy  at  the  source.  Ruzaikin  et  al 

show,  however,  that  the  L^  phase,  together  with  the  R^  phase,  can  be  drastically  j 

attenuated  by  structures  associated  with  the  Tibetan  Plateau.  The  result  is 

i 
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Earthquake*  Uaad  for  Hagnlcuda  Study 


Earthquake 

(HOAX) 

(Thie  Study) 

Dlf f trance 

♦tooth  Corolla* 

3.45 

3.59  i  0.19 

-0.14 

Hoy  19,  1971 

o  •  2 

o-6 

•SMs  Virginia 

4.5  1  0.4 

4.60  t  0.30 

-0.10 

lot.  20,  1999 

o  •  11 

n  -  13 

Tmumbm 

3.5 

3.96  t  0.17 

-0.46 

July  13,  1949 

o  •  4 

o*6 

••fir fiaia 

3.57  t  0.34 

3.51  t  0.17 

40.06 

May  31,  1944 

o-8 

n  ■  6 

••Hoot  Virginia 

3.62  l  0.1ft 

3.42  t  0.1ft 

♦0.20 

*rv.  24,  1994 

o-9 

n  -  4 

+* Alabama 

4.3  i  0.7 

3.65  t  0.7ft 

♦0.65 

fob.  11,  1964 

o  •  10 

a  -  6 

*  *  Included  la 

Ml 

■,  (MOM). 

**  U»ad  la  datazalnatlou  of  ft  cmly ;  other*  ueed  for  both  B  and  C 
determination* . 


AIT  versus  A  plou  for  three  southeastern  United  State*  earthquakes.  A  is  the  maximum 
amplitude  (microns)  of  the  L,  phase  on  the  short-period  vertical  seismogram  and  T  is  the  corresponding 
period  in  seconds.  A  is  the  epicentral  distance  in  degrees. 


Figure  13.  Event  magnitude  and  amplitude-distance  data  from  Bollinger  (1973). 
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Events  whic  i  Occurred  in  the  Southeastern  United  States  between  19G3  and  1975 


Event 

Date 

Latitude  mS 

Longitude 

*W 

Origin 
Time  UTC 
h  m 

No.  ol 
Stations 

Calculated 

Magnitude 

Reported  Magnitude 

i 

4  May  1963 

32.2 

79.7 

21 

01 

2 

3.3 

± 

0.02 

2 

b  Oct.  1963 

33.92 

82.51* 

06 

01 

3 

3.2 

± 

0.25 

3 

28  Oct.  1963 

36.7 

81.0 

22 

38 

2 

2.9 

± 

0.08 

4 

17  Feb.  1964 

34.7 

85.4* 

22 

47 

2 

3.3 

± 

0.21 

5 

18  Feb.  1964 

34.8 

85.5 

09 

31 

4 

4.2 

± 

0.10 

4.4  mi 

6 

7  Jii&r.  1964 

33.82 

82.5* 

13 

03 

3 

3.3 

+ 

0.15 

7 

13  Mar.  1964 

32.9 

83.4 

01 

20 

3 

3.9 

± 

0.15 

4.4  mi 

8 

20  Apr.  1964 

34.0 

81.1 

19 

05 

2 

3.5 

± 

0.04 

9 

12  Jun.  1965 

34.5 

83.8* 

10 

30 

2 

1.9 

0.05 

10 

20  Jun. 1965 

34.3 

87.1* 

12 

30 

3 

2.6 

± 

0.21 

11 

9  Sep.  1965 

33.5 

80.8 

14 

42 

3 

3.9 

± 

0.07 

12 

10  Sep.  1965 

33.5 

80.8 

02 

32 

3 

3.0 

± 

0.23 

13 

12  Sep.  1965 

34.8 

81.3 

18 

25 

4 

2.9 

0.27 

14 

8  Nov.  1965 

34.2 

82.8* 

12 

57 

3 

?..o 

± 

0.05 

15 

31  May  1966 

37.6 

78.0 

06 

19 

2 

4.1 

+ 

0.12 

3.1  mb 

16 

29  Jun.  '.967 

33.6 

90.9 

13 

57 

2 

3.6 

± 

0.07 

3.4  mb 

17 

23  Oct.  1967 

33.11 

80.71 

09 

04 

4 

3  9 

± 

0.21 

3.8  mb 

18 

10  Feb. 1968 

36.54 

89.90 

01 

34 

4 

4.3 

± 

0.08 

3.8  mb 

19 

8  Mar.  1968 

37.0 

80.5 

05 

38 

2 

4.1 

± 

0.18 

3.9  mb 

20 

22  Sep.  1968 

34.0 

81.5 

21 

41 

2 

3.5 

± 

0.07 

3-7  mb 

21 

9  May  1969 

33.92 

82.5* 

12 

13 

2 

3.1 

0.14 

22 

13  Jul.  1969 

36.4 

83.7 

21 

51 

4 

4.2 

0.10 

3.5  mi 

23 

20  Nov.  1969 

37.4 

81.0 

01 

00 

2 

4.7 

0.19 

4.8  ML,  CCS 

24 

11  Dec.  1969 

37.8 

77.4 

23 

44 

2 

3.4 

±1 

0.06 

4.3  mb 

25 

13  Dec.  1969 

35.08 

83.04 

10 

19 

4 

3.7 

0.15 

26 

30  Jui.  1970 

37.01 

82.25 

18 

48 

5 

3.3 

0.17 

3.8  mb 

27 

30  Jul.  1970 

37  .01 

82.25 

15 

15 

5 

3.7 

=fc 

0.11 

4.0  mb 

28 

31  Jul.  1970 

37.7 

03.4* 

00 

31 

2 

3.5 

± 

0.14 

29 

10  Sep.  1970 

36.1 

81.4 

01 

41 

4 

3.2 

± 

0.13 

30 

17  Nov.  1970 

35.9 

89.9 

02 

13 

3 

4.5 

± 

0.06 

3.6  mi 

31 

14  Mar.  1971 

33.1 

87.9 

17 

27 

2 

3.6 

± 

0.29 

3.9  Ml,  NOS 

32 

16  Mar.  1971 

33.1 

87.9 

02 

37 

3 

3.3 

± 

0.19 

33 

19  May  1971 

33.3 

80.56 

12 

54 

2 

3.9 

0  10 

3.4  mi 

34 

29  May  1971 

36.0 

82.0* 

21 

21 

2 

2.9 

± 

0.08 

35 

10  Jun.  1971 

34.7 

82.9* 

04 

19 

3 

2.8 

± 

0.19 

36 

13  Jun.  1971 

33.7 

86.6* 

01 

17 

2 

2.9 

± 

0.03 

37 

13  Jun.  1971 

36.0 

83.9 

02 

02 

2 

3.0 

± 

0.04 

38 

13  Jul.  1971 

34.7 

82.95 

09 

39 

2 

2.8 

± 

0.13 

39 

13  Jul.  1971 

34.7 

82.95 

10 

54 

2 

2.9 

± 

0.12 

40 

13  Jul.  1971 

34.7 

82.95 

11 

07 

2 

2.7 

±. 

3.02 

41 

13  Jul.  1971 

34.7 

82.95 

11 

42 

3 

4.4 

0.10 

42 

13  Jul.  1971 

34.7 

82.95 

11 

49 

3 

2.9 

± 

0.06 

43 

13  Jul.  1971 

34.7 

82.95 

15 

06 

2 

3.0 

± 

0.17 

44 

31  Jul.  1971 

33.7 

80.66 

20 

16 

4 

4.0 

± 

0.16 

45 

12  Sep.  1971 

38.1 

77.4 

00 

06 

2 

3.4 

± 

0.08 

46 

1  Oct.  1971 

35.8 

90.4 

18 

49 

2 

3.9 

± 

0.08 

4.11  Mia, 

47 

9  Oct.  1971 

35.9 

83.5 

16 

43 

2 

4.0 

± 

0.07 

3.4  mi 

48 

22  Oct.  1971 

36.0 

83.0 

21 

55 

3 

3.3 

dr 

0.11 

49 

1  Feb.  1972 

36.25 

90.9 

05 

42 

3 

4.0 

+ 

0.06 

4.1  mi 

3.9  M a,  SLM 

50 

3  Feb.  1972 

33.5 

80.4 

23 

11 

4 

4.3 

± 

0.16 

4.5  mb 

51 

7  Feb.  1872 

33.46 

80.58* 

02 

46 

4 

3.2 

± 

0.02 

Figure  14.  Event  list  and  amplitude-distance  data  from  Jones  and  Long  (1977). 
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Event 

Date 

Latitude  *N 

Longitude 

*W 

Origin 
Time  UTC 
h  m 

N'o.  of 
Stations 

Calculated 

Magnitude 

Reported  Magnitude 

52 

7  Feb.  1972 

33.46 

80.58’ 

02 

53 

4 

3.2 

± 

0.08 

53 

29  Mar.  1972 

36.2 

89.6 

20 

38 

2 

4.1 

+ 

0.05 

3.7  ms 

54 

5  Sep.  1972 

37.6 

77.7 

16 

00 

2 

3.6 

± 

0.04 

3.3  Mn„  BLA 

55 

25  May  1973 

33.9 

90.8 

14 

40 

2 

2.7 

± 

0.04 

56 

27  Oct.  1973 

28.7 

81.0 

06 

21 

2 

3.5 

± 

0.17 

57 

30  Oct.  1973 

35.7 

84.0 

22 

58 

4 

3.5 

lb 

0.13 

3.4  Mn„  SLM 

58 

30  Nov.  1973 

35.8 

84.0 

07 

48 

5 

4.5 

± 

0.16 

4.6  M, i(,  BLA 

59 

14  Dec.  1973 

35.7 

83.8* 

20 

58 

4 

3.1 

± 

0.08 

60 

19  Dec.  1973 

33.0 

80.3 

10 

16 

2 

3.0 

± 

0.12 

61 

8  Jan.  1974 

36.2 

89.39 

01 

12 

3 

4.4 

± 

0.04 

4.3  Mn„  SLM 

62 

12  Mar.  1974 

35.66 

89.79 

12 

30 

2 

3.3 

+ 

0.09 

4.1  nt 

3.2  Mn„  SLM 

63 

13  May  1974 

36.71 

89.39 

06 

52 

2 

4.1 

± 

0.01 

4 . 1  MbLg ,  SLM 

3 . 8  TTlb  a m 

64 

30  May  1974 

37.38 

80.42 

21 

28 

2 

3.7 

± 

0.22 

3.6  Mi i„  BLA 

65 

2  Aug.  1974 

33.87 

82.49 

08 

52 

4 

4.3 

± 

0.12 

4.9  MkL» i  GS 

4.3  rrtb 

4.8  MkL:  SLM 

66 

22  Nov.  1974 

32.9 

80.15 

05 

25 

3 

4.4 

± 

0.08 

4.7  mb 

67 

10  Dec.  1974 

31.3 

87.5 

06 

01 

4 

3.5 

=b 

0.02 

3.0  Mu,,  SLM 

68 

25  Dec.  1974 

35.78 

90.01 

13 

21 

2 

3.0 

± 

0.03 

3.0  ML 

69 

1  Mar.  1975 

33.55 

87.99 

11 

50 

2 

4.1 

dr 

0.26 

3.2  Mn„  SLM 

70 

1  Apr.  1975 

33.2 

83.2’ 

21 

09 

2 

3.9 

zb 

0.43 

71 

24  Jun.  1975 

33.72 

87.84 

11 

11 

2 

3.7 

rir 

0.03 

4.5  mt 

72 

29  Aug.  1975 

33.82 

86.6 

04 

22 

4 

4.4 

0.27 

4.4  MbL»t  SLM 

3.5  mb 

Figure  14(Cont.).  Event  list  and  amplitude-distance  data  from  Jones  and  Long 
(1977). 


50' 


that  earthquake  signals  which,  pass  through  such  structures  are  recorded  looking 
like  explosions.  For  example,  comparison  of  Figures  15a-d  from  Ruzaikin  et  al 
shows  this  effect.  Signals  1-11  are  recordings  similar  to  those  seen  in  the 
EUS  from  WUS  earthquakes.  Signals  13  through  15  are  similar  to  typical  explo¬ 
sion  records,  and  Signals  16  through  19,  although  they  are  from  earthquakes, 
are  more  explosion-like  than  those  from  any  explosion  we  have  observed  in 
the  United  States. 


Figures  from  Ruzaikin  et  al  show  the  character  of  received  signals  at 
stations  TLG,  NSB,  and  BDN.  Figure  15d  from  Ruzaikin  et  al  shows  some  of  the 
obstacles  to  efficient  propagation.  Figures  similar  to  these  will  be 
required  for  accurate  application  of  the  phase  as  a  discriminant. 

Bollinger  (1977)  has  measured  amplitudes  and  periods  for  the  list  of 
events  shown  in  Figure  16.  He  used  Nuttli's  (1970)  formula  to  compute  station 


magnitudes  and  then  plotted  the  normalized  amplitudes  as  a  function  of  distance. 
For  4°  to  30°  we  have  fitted  a  least-squares  slope  of  -2.1  as  indicated  in  the 


figure.  At  shorter  distances,  the  slope  seems  less,  but  again,  one  must  worry 
about  clipping  and  missed  maximum  cycles  at  such  distances  for  events  of 


magnitude  above  4.0  and  there  are  6  such  events  in  Bollinger's  data  set. 


Tables  I,  II,  and  III  give  a  summary  of  the  amplitude-distance  exponents 
implied  by  the  papers  reviewed  thus  far,  together  with  the  exponents  which  we 
have  determined  in  this  paper  from  Western  United  States  and  Eastern  United 
States  data.  In  general  the  results  found  in  this  study  are  compatible  with 
those  found  by  earlier  workers. 
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Map  of  Asia  showing  stations  used;  epicenters  of  earthquakes  for  which  seismograms 
are  shown  in  Figures  2  (Arabic  numerals),  3  (Roman  numerals),  and  1  (A  and  B);  and  areas 
of  high  mountains  (hatched  regions). 


Figure  15a.  L  propagation  characteristics  in  Asia  according  to  Ruzaikin  et 
a?  (1977). 
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Seismogram  of  selected  events  (Arabic  numerals  in  Figure  1)  recorded  at  a 
temporary  station  between  Talgar  (TLG )  and  Novosibirsk  (NSB) .  Events  1-9  are  a  profile 
extending  east  to  Baikal,  and  events  10-19  a  profile  to  the  south.  Recordings  are  from  a 
narrow  band  system  between  l.ls  and  2.3.  Scale  at  upper  right  shows  1  min. 

Figure  15b,  L  propagation  characteristics  in  Asia  according  to  Ruzaikin  et 
a?  (1977). 
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Map  shoving  interpretation  of  selected  paths  to  Soviet  stations:  (a)  triangles, 
stations;  (b)  circles,  epicenters;  (c)  dense  hatching,  apparent  sharp  boundaries  of  Lg 
propagation;  (d)  open  hatching,  less  sharp  boundaries;  (e)  dashed  hatching,  presumed 
southern  boundary;  (f)  solid  line,  clear  Lg;  (g)  dashed  line,  no  Lg;  (h)  dotted-dashed  line, 
veak  Lg;  and  (i)  heavy  dark  line,  selected  faults  on  the  Soviet  tectonic  map. 


Figure  15d.  L  propagation  characteristics  in  Asia  according  to  Ruzaikin  et 


A(km) 


Parameters  of  the  Southeastern  U.S.  Earthquakes  for  which  Lg  Attenuation  Data  were 

Obtained  in  this  Study 


Date 

State 

Origin  Time 
(UTO 

Epicenter 

(°N.  °W) 

Nt 

Feb.  18,  1964 

GA 

09:32:12 

34.8-85.5 

4.0  +  0.2 

8 

Nov.  25,  1964 

WV 

02:50:07 

37.4-81.7 

3.6  +  0.2 

13 

Apr.  26,  1965 

WV 

15:26:20 

37.3-81.6 

3.5  +  0.2 

7 

May  31,  1966 

VA 

06:19:02 

37.6-78.0 

3.6  +  0.2 

9 

Oct.  23,  1967 

SC 

09:04:10 

33.4-80.7 

3.4  +  0.2 

9 

Jul.  13,  1969 

TN 

21:51:09 

36.1-83.7 

4.2  +  0.2 

8 

Nov.  20,  1969 

WV 

01:00:10 

37.4-81.0 

4.6  +  0.3 

13 

Dec.  13,  1969 

NC 

10:19:34 

35.1-83.0 

3.4  +  0.1 

8 

Sep.  10,  1970 

NC 

01:41:10 

36.1-81.4 

3.1 

4 

Mar.  14,  1971 

AL 

17:27:52 

33.1-87.9 

3.6  +0.2 

7 

May  19,  1971 

SC 

12:54:03 

33.3-80.6 

3.7  +  0.2 

7 

Jul.  13,  1971 

SC 

11:42:26 

34.8-83.0 

3.8  +  0.3 

7 

Oct.  9,  1971 

TN 

16:43:34 

35.9-83.5 

3.7  +  0.2 

7 

Feb.  3,  1972 

SC 

23:11:08 

33.5-80.4 

4.5  +  0.1 

8 

Nov.  30,  1973 

TN 

07:48:41 

35.8-84.0 

4.6 

4 

Aug.  2,  1974 

GA 

08:52:10 

33.9-82.5  . 

4.1 

3 

Nov.  22,  1974 

SC 

05:25:56 

32.9-80.1 

4.3  +  0.2 

5 

*  Nuttli  (1973)  Lq  formula  used  to  determine  the  individual  stations's  m». 
t  Number  of  stations  averaged  to  obtain  the  listed  m»’s. 


Figure  16.  L  amplitude-distance  relations  from  Bollinger  (1977). 
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TABLE  I 


Decay  Rate  for  P 

n 


Author 

Range  (km) 

.  ~X 

x  in  r 

Comments 

Romney  (1959) 

200-800 

3.0 

Low  values  at  1000  km,  if  these 
points  neglected  -2.0  works  to 

2000  km.  Profile  to  Texas  and 
on  to  Maine.  Blanca  and  Logan  (WHS) 

Romney  et  al  (1962) 

200-1000 

3.0 

EUS  and  WUS  from  GNOME  (WUS) 

Ryall  and  Stuart 

(1963) 

200-800 

3.3 

WUS,  NTS  to  Colorado 

Nersesov  (1964) 

1.6 

Southern  border  USSR,  very  erratic 
from  Nersesov  Figure  6 

Evernden  (1967) 

2.0 

EUS  8.5  km/sec  refractor 

Evernden  (1967) 

200-1000 

3.04 

WUS  7.9  km/sec  refractor 

1000-1800 

2.0 

WUS  8.5  km/sec  refractor 

This  study 

200-2000 

2.5 

EUS 

200-1000 

3-4  ? 

WUS 
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TABLE  II 


Author 


Decay  Rate  for  P 

max 

Range  (km)  x  in  r  X 


Comments 


Romnev  et  al 
(1962) 

500-1500 

4 

GNOME 

p 

Ryall  and  Stuart 
(1963) 

200-1000 

3.5 

This 

study. 

fit 

to 

P  data 

Press  (1964) 

9 

3.2 

This 

study , 

fit 

for 

derived  Q 

Nersesov  (1964) 

-  200-1500 

2.4 

This 

study, 

fit 

to 

P  curves 
g 

Nersesov  (1969) 

200-2000 

2.3 

This 

study. 

fit 

to 

P  curves 

max 

Nersesov  (1964) 

200-2000 

3.1 

This  study, 
earthquake) 

fit 

to 

P  (Pribaikal 

max 

Evernden  (1967) 

~  200-1000  ? 

3.0 

P 

g 

This  study 

200-2000 

2.5 

EUS 

200-1500 

3.0 

wus 
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TABLE  III 


Decay  Rate  for  L  (R  ) 

g  g  max 


Author 

Range  (km) 

-x 

x  in  r 

Comments 

Richter  (1958) 

200-600 

3.0 

WUS 

Romney  (1959) 

200-2000 

3.0 

WUS,  all  3  components,  Blanca,  Logan 

Romney  et  al 
(1962) 

500-3000 

2.0 

EUS  GNOME 

Romney  et  al 
(1962) 

500-1500 

4.0 

WUS  GNOME 

Press  (1964) 

? 

2.2 

WUS,  this  study  fit  for  derived  Q 

Nersesov  (1964) 

200-2000 

2.3 

USSR 

Baker  (1970) 

250-500 

1.55 

WUS 

Baker  (1970) 

250-500 

3.46 

WUS 

Baker  (1970) 

200-1100 

2.71 

WUS 

Baker  (1970) 

200-2000 

2.30 

WUS 

Nuttli  (1973) 

50-400 

0.9 

EUS 

Nuttli  (1973) 

400-3000 

1.66 

EUS  (slope  assumed  and 
not  rejected) 

Bollinger  (1973) 

300-3000 

2.0 

EUS,  cannot  reject  this  slope 

Street  (1976) 

400-3000 

-  2.5 

EUS,  fit  by  eye,  this  study.  St. 
Lawrence  earthquake 

Jones  and  Long 
(1977) 

400-900 

2.0 

EUS,  cannot  reject  this  slope 

Bollinger  (1977) 

150-3000 

2.0 

EUS,  cannot  reject  this  slope 

This  study 

250-2000 

2.0 

EUS 

This  study 

200-1500 

3.0 

WUS 
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I ' i scr  Lmination 

Will is,  DeNoyer  and  Wilson  (1963)  and  Willis  (1963)  showed  that,  although 
there  is  some  overlap,  the  ratio  of  the  maximum  of  the  "shear-type"  waves  to 
the  maximum  of  the  first  10  seconds  or  so  of  the  compressional  motion  is 
greater  for  earthquakes  than  for  explosions.  This  is  as  recorded  at  distances 
less  than  10°  and  for  frequencies  from  0.5  to  10  Hz.  Distance  corrected  net¬ 
work  averages  are  not  discussed  in  these  studies;  only  population  means  or 
single-station-event  measurements  as  a  function  of  distance  or  frequency. 

Booker  and  Mitronovas  (1964)  used  stations  in  the  general  range  of  15° 
to  define  average  energy  ratios  of  different  phases  (velocity  windows).  Fair 
discrimination  was  achieved  using  ratios  of  shear  to  compressional  phases. 

This  is  the  only  study  in  the  literature  to  use  network-averaged  short-period 
discriminants . 

Lambert  and  Becker  (1975)  used  vertical  records  of  regional  events  re¬ 
corded  at  NORSAR  to  show  that  there  is  some  suggestion  that  the  shear  to 
compressional  ratio  is  greater  for  earthquakes  than  for  explosions. 

Bakun  and  Johnson  (1970),  using  a  rather  small  sample  of  natural  earth¬ 
quakes,  found  that  at  a  regional  station  natural  earthquakes  have  less  high 
frequency  in  the  phase  than  do  NTS  explosions.  However,  aftershocks  of  the 
explosions  were  similar  to  the  earthquakes.  Other  workers,  Peppin  and  McEvilly 
(1974),  Murphy  (1976)  found  no  spectral  discrimination.  Peppin  (1977  personal 
communication)  has  said  that  use  of  Bakun  and  Johnson’s  discriminant  failed  on 
the  "natural"  Massachusetts  Mountain  earthquakes  at  NTS  south  of  Yucca  Flats. 

Pasechnik  (1970)  showed  travel  times  for  P  and  S  phases  out  to  200  km  and 
asserted  that  S/P  amplitude  ratios  are  less  for  explosions  than  for  earthquakes. 
He  also  plotted  P  wave  periods  for  earthquakes  and  explosions  as  a  function  of 
distance  in  the  range  250  to  800  km.  The  explosion  periods  are  strikingly 
shorter.  This  may  be  due  to  earthquakes  being  located  in  tectonic  regions 
while  explosions  are  located  in  stable  regions. 
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DISTANCE  AMPLITUDE 


Tables  IV  and  V  are  lists  of  the  earthquakes  and  explosions  analyzed  in 
this  study  for  the  purpose  of  determining  accurate  amplitude-distance  relation¬ 
ships  and  to  determine  the  ability  to  use  regional  phases  for  discrimination. 

One  of  the  most  illustrative  events  is  the  Hebgen  Lake  earthquake  of  1964.  This 
event  is  almost  exactly  on  the  border  between  the  Western  United  States  (WUS) 
and  Eastern  United  States  (EUS)  "absorption  regions"  as  characterized  by 
Der,  Massd  and  Gurski  (1975).  See,  e.g.,  their  Figures  12  and  13. 

Thus  waves  from  this  event  offer  an  excellent  opportunity  to  see  the 
difference  in  attenuation  between  the  EUS  and  WUS.  In  fact,  in  Figure  17 
we  see  that  a  decay  rate  of  -3  is  appropriate  for  L^  in  the  WUS  and  a  value 
of  -2  in  the  EUS.  We  also  see  very  vividly  how  the  first  arrival  on  the  south¬ 
east  profile,  at  GVTX,  has  a  low  amplitude  which  may  be  understood  by  virtue 
of  the  fact  that  this  path  is  along  the  EUS-WUS  boundary  as  delineated  by 
Der  et  al  (1975).  The  large  value  at  NPNT  may  be  partially  understood  in  that 
the  measured  period  at  this  station  was  2  seconds,  leading  to  a  large  period 
correction  for  the  amplitude.  Note  that  the  EUS  decay  seems  to  be  similar  to 
that  in  the  USSR. 

In  Figure  18  we  see  a  similar  plot  for  P  and  see  that  the  difference 

max 

between  EUS  and  WUS  is  not  as  large  as  for  L  ;  the  data  are  consistent  with 

8 

-3  for  WUS  and  -2.5  for  EUS.  Note  that  the  P  decay  for  the  EUS  seems  to 

max 

be  similar  to  that  in  the  USSR. 

In  Figure  19  we  see  the  amplitude  as  a  function  of  distance  for  both  L  and 

8 

^max  ^or  t*ie  event  SALMON  and  the  closely  matched  18  February  1964  Alabama 
earthquake.  Signals  from  SALMON  passed  almost  directly  through  the  epicenter 
of  this  earthquake,  only  5°  to  the  Northeast  on  the  way  to  stations  BLWV,  BRPA, 
DHNY,  LSNH,  and  HNME.  A  particularly  well-matched  set  of  stations  are  BRPA, 
approximately  8°  from  the  earthquake,  and  BLWV  approximately  10°  from  SALMON. 

We  shall  make  extensive  use  of  the  recordings  of  SALMON  at  BLWV  and  18  February 
at  BRPA  in  comparisons  of  the  two  events. 

As  an  indication  of  the  discrimination  capabilities  of  the  two  phases  we 

see  in  Figure  19  that  although  the  Lg  amplitudes  are  of  equal  level  for  the  two 

events  the  P  amplitude  for  SALMON  is  on  the  average  0.3  magnitude  units  (mu) 

above  that  of  18  February.  We  see  also  from  this  figure  that  the  slopes  of  2.0 

and  2.5  for  L  and  P  seem  appropriate, 
g  max 
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TABLE  VII 


A/'T  (m^i/sec) 


La  HEBGEN  LAKE 


A/T  (m/i/sec) 


Since  the  event  SALMON  is  ol  special  interest  as  the  only  KUS  exp  I  os i  mi , 
we  shall  include  in  the  body  of  this  report,  the  P  and  Lg  plots  of  amplitude 
data  for  all  profiles  out  of  the  epicenter.  The  P  data,  separated  into  first 
arrival  P  and  Pmax  readings,  is  found  in  Figure  20.  A  least-squares  slope  of 
-2.54  has  been  fitted  through  the  Pmax  data  from  the  Northeast  and  North 
profiles.  We  see  that,  in  general,  data  at  greater  than  9°  distance  along 
the  Western  profile  falls  well  below  this  line.  In  Figure  21  we  have  the 
SALMON  Lg  data.  Analysis  of  this  phase  for  this  event  must  keep  in  mind  that 
in  many  cases  a  Mexican  event  is  mixed  in,  thus  raising  the  amplitude.  A 
least-squares  fit  to  the  uncontaminated  Northeast  and  North  profile  data  yields 
a  slope  of  -2.03.  Note  the  large  reduction  of  amplitude  between  WUS  and  EUS 
stations.  This  difference  is  greater  than  that  seen  for  the  P  phases  and 
suggests  that  complications  in  structure  affect  L  more  than  P.  This  would 

O 

then  suggest  that  complications  in  structure  are  apt  to  result  in  earthquakes 
being  identified  as  explosions  unless  a  careful  regionalization  is  undertaken. 

The  distance-amplitude  relations  of  the  last  event  to  be  examined  in  the 
main  body  of  the  report  is  that  of  HARDHAT.  This  shot  was  in  granite  in  the 
WUS  and  was  comprehensively  recorded.  We  see  in  Figure  22  that  the  P  data  is 
extremely  erratic.  However,  the  Pmax  data  is  more  stable  and  a  regression  line 
slope  of  -3,37  has  been  fitted.  This  would  be  appropriate  for  the  WUS.  In 
Figure  23  we  see  the  data  for  the  HARDHAT  Lg  phase.  Here  a  slope  of  -3.13 
has  been  determined  by  least  squares.  Stations  HBOK  and  WMO  in  the  EUS  have 
unusually  high  amplitudes,  as  might  be  expected  if  the  decay  rate  changed  to  a 
reduced  value  when  the  waves  emerged  from  the  WUS. 

In  Appendix  I  we  have  reproduced  distance-amplitude  figures  similar  to  those 

just  discussed  for  events  of  special  interest,  GNOME,  four  EUS  earthquakes,  SHOAL, 

the  earthquake  FALLON,  and  BILBY.  Consideration  of  the  least  squares  slopes 

determined  as  discussed  above  together  with  the  historical  results  determined  by 

other  researchers  as  tabulated  in  Tables  I,  II  and  III  suggests  the  hypothesized 

slopes  of  2.5,  2.5  and  2.0  in  the  EUS  for  P  ,  P  ,  and  L  respectively.  For  the 
K  ’  n’  max  g 

WUS  the  hypothesized  slopes  are  3.5,  3.0  and  3.0.  In  Table  VI  we  see  that  in 
general  these  slopes  cannot  be  rejected  for  the  data  which  have  been  plotted; 
and  so,  we  adopt  these  values  for  our  discrimination  analysis. 
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Figure  20.  Amplitude-distance  relations  for  the  P  and  P  (phase  velocity 
greater  than  5  km/sec)  phase  from  SALMON.  Amieast-squares  slope 
of  -2.54  has  been  fitted  to  the  PMX  data  of  the  Northeast  and 
North  profiles.  Note  that  most  of  the  data  along  the  Western 
profile  fall  well  below  the  line. 


Northeast  •  EUAL.  CPO.  BLWV.  DHNY  LSNH.  HNME 

JELA.  GVTX,  WMO,  FOTX  RTNM,  DRCO. 
KNUT.  EKNV,  HL2IO 

North  4  VOIO.  WFMN.  GPMN,  RYND 
TSND,  RKON 

V 

'  Lg  (Rg)  max  for  Salmon 


*^N)(M)  M  (Mixed  with  Mexican  Earthguake) 
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Figure  21.  Amplitude-distance  relation  for  L  (R  )  for  SALMON.  A  least- 

squares  slope  of  -2.03  has  been  fftt§dmto  the  data  along  the 
Northeast  and  North  profiles.  Note  that  most  of  the  data  along 
the  Western  profile  fall  well  below  the  line.  Some  of  the  data 
has  been  contaminated  by  signals  from  an  earthquake  in  Mexico 
which  occurred  this  day.  These  measurements  are  marked  by  (M) . 
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Figure  22.  Amplitude-distance  relations  for  the  P  and  P  (phase  velocity 
greater  than  5  km/sec)  phase  from  HARdRaT.  Amfeast-squares  slope 
of  -3.37  has  been  fitted  to  the  P  data.  Note  that  the  P  data 
is  much  more  erratic.  The  notati§nX(c)  indicates  that  the  9ata 
were  clipped  and  that  the  true  value,  therefore,  lies  above  the 
piotted  point,  (N7)  indicates  that  the  measurement  may  be  noise. 
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Amplitude-distance  relation  of  L  (R  )  for  HARDHAT.  A  least- 
squares  slope  of  -3.13  has  been  ?it£e3a£o  the  data. 


TABLE  VI 


Summary  of  amplitude-distance  slopes 
EUS 


P 

n 

P 

max 

L 

g 

SALMON 

>  2.5 

2.54  -  2.5 

2.03  + 

2.0 

New  Madrid 

2.5  OK 

2.5  OK 

2.0  OK 

Massachusetts  Coast 

2.5  OK 

2.5  OK 

2.0  OK 

Southern  Illinois 

2.5  OK 

2.5  OK 

2.0  OK 

New  Hampshire 

2.5  OK 

2.5  OK 

2.0  OK 

GNOME 

2.5  or  > 

2.5  OK 

2.5  OK 

2.0  OK 

WUS 

Pn 

pmax 

Lg 

HARDHAT 

3.5  ? 

3.37  -+  3.0 

3.13  -*• 

3.0 

SHOAL 

4.0 

2.77,  3.0  OK 

3.0  OK 

FALLON  emergent , 

erratic 

2.0,  3.0  poor 
fit 

3.0  OK 

BILBY 

3.5  OK 

2.0,  3.0  OK 

3.0  OK 

GNOME 

3.0  OK 

3.0  OK  or  » 
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very  small 
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DISCRIMINATION  -  P  :L 
max  g 

In  the  discrimination  analysis  we  measured  periods  and  amplitudes  for 

P  and  L  for  each  of  the  events  in  Tables  IV  and  V.  Measurements  were 
max  g 

made  at  the  number  of  stations  indicated  in  the  tables.  Care  was  taken  that 
the  stations  were  as  well  distributed  as  possible  and  that  the  magnitudes  of 
the  TFO  extended  array  in  Arizona  were  averaged  together  to  give  an  array 
magnitude  which  was  then  averaged  in  with  the  other  stations.  Magnitudes  were 
computed  only  from  paths  for  which  at  least  the  first  2/3  of  the  path  was  in 
the  same  province  (EUS  or  WUS)  as  the  event.  The  events  Hebgen  Lake  and 
GNOME  were  considered  to  be  on  the  EUS-WUS  boundary  so  that  all  paths  were 
in  one  province  or  the  other. 

Magnitudes  were  calculated  to  be  the  measured  amplitude  corrected  to  1000  km 
by  the  appropriate  formula  and  averaged  over  the  network.  The  results  are  seen 
in  Figure  24  where  the  numbers  are  keyed  to  Table  IV.  First,  we  see  that  the 
EUS  earthquakes  are  well  separated  from  SALMON,  the  only  EUS  explosion;  and 
from  the  GNOME  event  when  its  magnitude  is  computed  only  on  the  EUS  paths.  The 
separation  between  the  population  lines  of  slope  1.0  is  0.6  magnitude  units. 

For  WUS  events  the  situation  is  more  complicated.  The  discrimination  is 

not,  on  the  face  of  it,  as  good  as  in  the  EUS.  If  we  simply  take  the  difference 

between  the  population  means,  then  we  see  a  difference  of  approximately  0.45 

magnitude  units,  slightly  less  than  the  0.6  in  the  EUS.  However,  while  the 

earthquake  population  shows  only  slightly  greater  scatter  than  in  the  EUS,  the 

scatter  in  the  explosion  population  is  quite  large.  However,  the  principal 

events  which  contribute  to  this  scatter  are  BILBY,  MONERO,  YUBA  and  PLATTE. 

BILBY  is  so  large  that  its  low  corner  frequency  may  have  caused  the  P  ^  signal 

to  lose  amplitude  relative  to  L  because  the  L  always  is  of  low  frequency  due  to 

6  8 

modal  cut-off. 

YUBA  and  PLATTE  were  selected  because  they  were  extremely  small  events 
and  MONERO  because  it  was  the  closest  NTS  event  to  the  Massachusetts  Mountain 
Earthquake.  YUBA  and  PLATTE  are  at  a  depth  of  only  about  200  meters  and  were 
detonated  in  a  tunnel  at  Rainier  Mesa.  These  unusual  conditions  may  be 
responsible  for  a  low  corner  frequency  for  P  ,  leading  to  a  low  amplitude. 
Evidence  for  this  is  that  Pn  at  KNUT  has  a  period  of  0.2,  0.6  and  0.4  seconds  for 
BUTEO,  YUBA,  and  PLATTE  respectively.  BILBY  has  a  period  of  0.6  seconds;  KNUT 
was  not  up  for  MONERO.  BUTEO  is  probably  more  characteristic  of  a  true  test  ban 
monitoring  situation  in  that  it  was  overburied,  a  procedure  necessary  to  avoid 
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ilc to  it  ion  of  cratering  phenomenon.  If  we  do  take  the  remaining  WUS  explosions 
to  be  the  population  of  interest  the  explosion  population  is  almost  as  tight 
as  the  earthquake  population  and  the  separation  is  about  0.55  magnitude  units, 
almost  identical  to  the  0.6  in  the  WTJS. 

Examining  the  events  Hebgen  Lake  and  GNOME  in  the  EUS  and  WUS  shows  that 

for  L^  the  WUS  amplitude  is  about  0.3  to  0.55  magnitude  units  less  at  1000  km 

than  is  the  EUS  amplitude.  Let  us  take  0.42  magnitude  units  as  the  mean  dif- 

-2  -3 

ference.  Then,  assuming  r  in  the  EUS  and  r  in  the  WUS  the  waves  would  be 

of  equal  amplitude  at  380  km.  This  can  be  seen  to  be  a  reasonable  conclusion 

by  examining  Figures  17  and  Al-lb.  This  then  further  suggests  that  somewhere 

around  380  km  the  amplitude-distance  curves  of  EUS  and  WUS  should  begin  to  merge 

and  this  may  offer  a  start  toward  a  unified  magnitude  scale  for  the  entire  United 

States.  In  this  connection  we  should  note  that  Nuttli  (1973)  showed  that  y  =  0.6 

adequately  fitted  Richter's  (1958)  magnitude  relation  between  40  km  and  700  km. 

-3 

It  is  also  true  that  a  slope  of  r  closely  fits  Richter’s  curve  between  200  and 
600  km. 

Discrimination-Spectra  and  Transverse/Radial 

Figures  25  and  26  show  respectively  the  Alabama  earthquake  of  18  February 
as  recorded  at  BRPA,  a  distance  of  823  kilometers,  and  SALMON  as  recorded  at 
station  BLWV  at  a  distance  of  1058  km.  These  two  stations  are  indicated  by 
Der,  McElfresh  and  Mrazek  (1979)  to  have  the  same  crustal  structure  so  that 
any  differences  between  the  recordings  are  probably  not  due  to  differences  in 
the  crustal  structure. 

We  see  that  the  transverse  to  radial  maximum  amplitude  ratios  before  and 

after  5  km/sec  are  the  same  for  the  earthquake  and  explosion,  and  that  in  the 

time  domain  there  is  no  noticeable  difference  in  the  frequency  content  between 

components,  or  between  events  for  common  phases.  One  distinctive  difference  which 

can  be  seen,  however,  is  that,  especially  as  seen  on  the  transverse  component, 

the  earthquake  Lg  begins  abruptly  and  grows  quickly  to  a  maximum  amplitude. 

In  contrast,  the  SALMON  L  shows  no  clear  arrival  and  increases  slowly  to  its 

g 

maximum.  These  differences  might  be  explained  by  means  of  decreased  SALMON 
amplitude  of  shear  waves  which  are  needed  for  a  dramatic  emergence  of  the 
L  from  the  P  and  S  coda,  or  by  the  (presumably)  deeper  focus  of  the 
earthquake  which  excites  a  higher  proportion  of  higher  modes  which  have  faster 
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group  velocities.  Other  possible  influences  are  the  difference  in  the  crust 
at  the  source  and  the  fact  that  the  earthquake  signal  does  not  propagate 
through  the  thickly  layered  sediments  which  lie  between  the  SALMON  and 
earthquake  epicenters.  The  above  observations,  are,  therefore,  mainly  of  use 
as  guidelines  for  further  research. 

Another  point  which  is  easily  seen  in  Figure  25  on  the  transverse  com¬ 
ponent  is  that  there  is  a  sharp  decrease  in  frequency  with  the  arrival  of  the 
phase.  This  figure  also  shows,  however,  that  it  will  be  difficult  to  obtain 
a  "true"  L^  spectrum  at  high  frequencies  because  immediately  preceding  its 
arrival  there  is  coda  signal,  and  this,  presumably,  continues  "under"  the  L 

8 

signal.  Thus,  even  if  there  is  a  modal  cutoff  for  L  at  high  frequencies,  there 

will  be  high  frequency  energy  in  the  window  which  is  greater  than  the  ambient 

noise  before  the  arrival  of  P  .  Great  care  must  therefore  be  taken  in  inter- 

n 

preting  this  energy  in  terms  of  the  source  spectrum  via  an  L^  normal  mode  theory 
which  does  not  take  account  of  scattering  from  earlier  phases.  We  shall  return 
to  this  question  by  using  spectral  analysis. 

The  major  point  to  make  with  respect  to  discrimination  is,  of  course,  that 

the  L  / P  ratio  is  greater  for  the  earthquake  than  for  SALMON.  This  is  easily 

g  max  J 

seen  on  all  three  components,  and  reference  can  be  made  to  Figure  19  to  see  that 
the  Lg/Pmax  ratio  is  about  a  factor  of  2  greater  for  the  earthquake  than  for 
SALMON.  This  is  less  than  the  average  separation  of  ~  0.6  magnitude  units  seen 
in  Figure  24. 


Finally,  P^  is  much  more  emergent  for  the  earthquake  than  for  the  explosion. 

The  P  /P  ratio  is  about  1:1  for  SALMON  and  1:4  for  the  earthquake.  This 
n  max 

point,  however,  has  not  been  studied  in  detail  and  may  be  closely  related  to 
crustal  structure  so  we  cannot  assert  that  the  Lg/Pn  ratio  would  be  a  reliable 

discriminant. 


Proceeding  to  the  spectral  details,  we  see  in  Figure  27a, b,c  the  Z,  R,  and  T 
SALMON  P-wave  signal  spectra  and  noise  at  BLWV,  and  in  Figure  28  the  corresponding 
signal  for  the  18  February  1964  earthquake  at  BRPA.  The  spectra  have  not  been 
corrected  for  the  LRSM  short-period  instrument  response.  In  the  time  domain, 
as  remarked  before,  the  signals  seem  to  have  similar  frequency  content  and  by 
overlaying  the  two  log  spectral  plots  for  vertical  components  one  can  see  that 
they  are  identical  to  within  experimental  error  as  indicated  by  the  dashed  line 
in  Figure  28.  Der  and  McElfresh  (1976)  have  shown  that  the  P  spectrum  of 
SALMON  at  BLWV  and  BRPA  is  almost  perfectly  matched  by  the  SALMON  reduced  dis¬ 
placement  potential  without  absorption.  Thus  the  spectrum  in  Figure  27a  is  the 
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Figure  27a.  SALMON  P  waves  and  noise,  signal  and  spectra  at  BLWV,  ver 
component. 
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SALMON  source  spectrum  as  seen  through  the  LRSM  instrument  response,  and  thus 
the  18  February  earthquake  has  a  source  spectrum  identical  to  that  of  SALMON 
between  1  and  10  Hz.  The  same  conclusion  is  reached  when  the  signals  are  com¬ 
pared  at  other  stations.  Comparison  of  the  vertical  and  transverse  components 
spectra  in  Figure  27c  shows  that  the  transverse  P  is  of  slightly  higher  fre¬ 
quency  than  the  vertical. 

The  similarity  of  the  earthquake  and  explosion  spectra  suggests  that 
spectral  discrimination  between  earthquakes  and  explosions  may  be  very  difficult, 
if  not  impossible.  It  may  be  objected  that  this  is  only  one  pair  of  events; 
however,  most  studies  of  spectral  differences  between  earthquakes  and  explosions 
use  explosions  not  in  tectonic  regions  so  that  differences  in  propagation  may 
explain  any  differences  seen.  Also,  those  spectral  differences  which  do  exist 
may  be  due  to  the  cancellation  of  explosion  energy  at  frequencies  well  below 
1  Hz  by  interference  between  pP  and  P.  Since  for  Figure  28  S/N  <  1  for  f  <  1  Hz, 
this  frequency  range  may  possibly  still  be  exploited  for  discrimination;  however, 
it  is  clear  that  in  practice  it  will  be  necessary  to  develop  methods  for  decreasing 
the  low-frequency  noise. 


If  the  P-wave  and  other  source  spectra  of  the  earthquake  and  explosion 

are  similar,  then  one  might  expect  that  the  L  spectral  shapes  would  be  similar 

8 

to  each  other  since  the  propagation  effects  and  modal  cutoffs  are  the  same. 
Comparing  Figures  29a  and  30  (see  overlay  on  Figure  30)  we  see  that  this  is, 
in  fact,  the  case.  By  comparing  Figure  29a  to  Figure  27a  or  30  to  28  we  see 
that  the  spectrum  is  shifted  to  lower  frequencies  as  compared  to  the  P  wave 
spectrum.  Since  the  P-wave  spectrum  is,  as  we  have  seen,  equal  to  the  source 
spectrum,  it  follows  that  the  L^  spectrum  is  not  the  source  spectrum.  This 
might  be  expected  because  of  the  existence  of  modal  cutoffs.  It  follows  that 
determinations  of  corner  frequencies  and  related  source  parameters  from  L  , 
see  e.g.  Street  (1976)  are  liable  to  lead  to  errors. 


Figure  29d  as  compared  to  Figure  27a  shows  that  the  P^  spectrum  is  almost 

the  same  as  P  .  P  has  somewhat  more  energy  around  2  Hz  and  is  therefore  some- 
n  g 

what  lower  in  overall  frequency  content. 


The  true  L^  spectrum  is  probably  even  more  low-frequency  than  is  immed¬ 
iately  apparent  by  comparison  of  Figure  27a  and  29a.  In  Figure  29a  we  have  also 

plotted  the  P  spectrum  seen  in  Figure  29d.  We  see  that  it  merges  into  the  L 
8  8 
spectrum  near  2.4  Hz.  Thus  the  energy  in  the  L^  spectrum  above  2.4  Hz  is 

probably  scattered  coda  energy  and  the  true  L  spectrum  may  be  plunging  very 
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Figure  30. 


18  February  L  waves  and  noise,  signal  and  spectra  at  BRPA 
note  superposition  of  SALMON  L  spectrum  from  BLWV. 


deeply  above  2.5  Hz.  It  will  be  necessary  to  keep  this  fact  in  mind  when 
comparing  observed  and  theoretical  spectra. 


The  superposition  of  the  vertical  spectrum  on  the  transverse  in  Figure  29c 

shows  that  the  transverse  is  slightly  lower  in  frequency  in  contrast  to  P  and 

n 

P  in  which  the  transverse  is  of  higher  frequency.  Comparison  of  Figures  29d,e,f 
shows  that  above  5  Hz  the  transverse  and  to  some  degree  the  radial  components  of 
P  have  somewhat  more  energy  than  the  vertical  . 


The  contamination  of  the  spectra  by  coda  energy  may  have  made  it  possible 

for  Street  (1976)  and  other  workers  to  determine  fairly  accurate  corner  frequency 

and  high-frequency  asymptotes  from  L  spectra  even  though  the  procedure  would  not 

8 

be  justifiable  from  a  strict  theoretical  point  of  view. 


Figures  31  to  34  enable  us  to  reach  conclusions  identical  to  those  above 
by  using  station  EUAL  which  was  242  km  from  SALMON,  and  311  km  from  the  18 
February  earthquake. 

Figures  35  to  38  are  included  for  completeness.  They  give  SALMON  as  seen 
at  BRPA,  and  the  18  February  earthquake  as  seen  at  BLWV.  Thus  the  events  may 
be  compared  at  common  stations  but  different  distances,  instead  of  at  more  equal 
distances  and  different  stations,  as  we  did  in  the  analysis  of  Figures  27  to  30; 
the  same  conclusions  are  reached. 

Finally,  SALMON  signals  of  special  importance  are  also  included:  the  P  and 
Lg  waves  as  seen  at  station  JELA,  a  distance  of  243  km,  (39,  40);  and  at  stations 
HNME  (41)  and  RKON  (41,  42).  These  stations  are  constantly  used  as  a  basis  of 
comparison. 

In  particular,  note  in  Figure  41  that  the  spectrum  at  HNME  has  less  high 
frequency  than  at  RKON,  even  though  they  are  at  comparable  distances  of  22.46° 
and  19.9°,  respectively.  The  frequency  content  of  the  signals  seems  to  be 
stationary  in  the  time  plots,  suggesting  that  triplications  are  not  altering  the 
general  trends  of  the  spectra.  Thus,  there  is  apparently  more  absorption  under 
HNME  than  under  RKON.  The  spectral  ratio  between  1  and  2  Hz  in  Figure  41 
suggests  a  At*  value  of  0.4  sec. 


Results  at  NOR5AR 

The  relative  frequency  content  of  P  and  L^  at  regional  distances  which  we 
have  just  noted  for  SALMON  and  for  the  February  18  earthquake  may  also  be 
observed  for  Soviet  explosions  as  observed  at  NORSAR.  In  Table  VII  is  a  list  of 
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events  which  we  shall  analyze.  The  first  two  events  in  that  list  are  the 

closest  shots  to  NORSAR  and  their  signals  and  spectra  are  given  in  Figures 

44  to  48.  Comparing  Figures  44  and  45  for  the  event  at  11.92°  we  see 

that  the  spectrum  is  shifted  to  lower  frequencies  as  compared  to  the  P 

spectrum,  in  agreement  with  what  we  saw  in  comparing  Figures  1.27a,  29a), 

(32,  33),  (37,  38)  and  (42,  43)  for  SALMON.  The  same  shift  is  evident  for 

a  Soviet  explosion  at  the  distance  of  15.69°  as  seen  in  Figures  46  and  47. 

(Figure  45  is  from  the  LP  system  at  NORSAR,  and  by  comparison  with  Figure  44 

we  see  that  it  is  far  easier  to  detect  L  at  2  Hz  than  LR  at  0.05  Hz, 

g 

approximately  26  dB  or  1.3  magnitude  units  easier.) 
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Figure  33.  SALMON  L  waves  and  noise,  signals  and  spectra  at  EUAL  a  distance 
of  242  kfi. 
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Figure  36. 


SALMON  L  wave  signals  and  noise,  waveforms  and  spectra  at  BRPA. 
Note  overlay  of  spectrum  of  coda  from  time  just  before  L  arrival 
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Figure  43.  P  wave  signal  and  noise,  waveform  and  spectra  from  Soviet 
explosion  on  72/09/04  at  a  distance  of  11.92°  as  recorded 
at  C3  element  of  NORSAR.  See  Table  VI  for  details. 


SIGNAL  NOISE 


NR33B2LG 
20002  3 

23  9  my 


2.5  sec 


(■ 


100.0 

(m|i2/Hz) 


- 10 

10  Hz 


Figure  44.  Lg  wave  signal  and  noise,  waveform  and  spectra  from  Soviet 
explosion  72/09/04  at  a  distance  of  11.92°  as  recorded  at 
C3  element  of  NORSAR.  See  Table  VI  for  details. 
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Figure  45. 


0.5  Hz 

Long  period  LR  wave  signal  and  noise,  waveform  and  spectra  from 
Soviet  explosion  72/09/04  at  a  distance  of  11.92°  as  recorded 
at  C3  element  of  NORSAR.  See  Table  VI  for  details. 
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COMMENTS  ON  DETECTION 

Inspection  of  Figures  (43,  44)  and  (46,  47)  reveals,  as  mentioned  before, 

that  the  P  waves  have  more  high  frequency  than  the  L  waves.  Considering 

S 

that  the  P  wave  amplitude  for  explosions  for  the  15.69°  event  is  larger  than 
the  L  ,  it  would  seem  that  detection  of  explosions  (not  necessarily  earth- 

s 

quakes)  will,  in  general,  hinge  on  the  P  waves.  It  is  possible  with  the 
spectra  of  the  signals  in  Table  VII,  as  presented  in  Appendix  II,  to  give  a 

rough  estimate  of  the  detection  threshold  of  a  single  instrument  at  the  site 
of  the  C3  subarray.  Let  us  assume  that  an  analyst  will  detect  the  P  wave  when 
the  signal  spectrum  is  equal  to  the  noise  spectrum  at  some  frequency.  This 
is  reasonable  for  two  reasons.  First,  the  signal  to  noise  ratio  is  typically 
flat  over  a  bandwidth  of  nearly  2  Hz.  This  is  wider  than  is  typically  the  case 
for  teleseisms  and  makes  the  "frequency  change"  more  dramatic.  Secondly,  these 
regional  signals  are  more  extended  in  duration  than  the  teleseismic  signals, 
thus  giving  the  analyst  a  larger  window  in  time  over  which  to  establish  that 
the  character  of  the  trace  has  permanently  changed  from  that  before  the  arrival 
of  the  signal.  (In  theoretical  signal  detection  terms,  these  remarks  are  equi¬ 
valent  to  the  statement  that  regional  P  signals  have  a  larger  BT  or  bandwidth¬ 
time  product  than  teleseismic  signals,  so  that  they  may  be  detected  at  a  lower 
threshold.) 

The  required  calculations  are  outlined  in  Table  VIII.  There  we  note  the 

S/N  in  dB  at  the  frequency,  f^  ,  at  which  the  S/N  is  a  maximum.  Assuming  that 

the  S/N  is  proportional  to  magnitude,  we  have  then  used  the  ISC  m^  values  to 

correct  these  dB  values  to  those  appropriate  for  a  magnitude  4.5  event.  These 

are  then  further  corrected  to  what  they  would  be  at  a  distance  of  10°  assuming 

-2.5 

amplitude  proportional  to  r  '  as  we  found  in  the  EUS  and  by  consideration  of 

Nersesov  and  Rautian's  (1964)  results.  These  dB  values  are  then  averaged  to 

-2.5 

give  a  value  of  43  dB  S/N  at  10°  for  an  m^  =  4.5  explosion.  Then,  using  r 
for  extrapolation,  we  find  m^  detection  thresholds  of  0.9,  2.4,  3.3,  and  3.9  m^ 
at  5°,  10°,  15°  and  20°  respectively.  These  are  probably  low  estimates  to  some 
degree,  because  such  small  events  would  have  more  high-frequency  energy  than 
the  events  considered  as  data.  We  have  not  considered  whether  this  is  the  50% 
or  90%  threshold  so  there  is,  perhaps,  a  residual  0.3  unit  uncertainty. 
However,  since  we  think  analysts  would  almost  certainly  detect  a  properly 
filtered  signal  under  these  conditions,  our  best  estimate  is  that  it  is  a 
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90%  threshold.  It  must  also  be  stated  that  this  is  a  shield  to  shield  result. 
Events  in  tectonic  regions  would  certainly  have  less  high-frequency  and  be  harder 
to  detect. 

It  is  worth  noting  that  although  in  Table  VIII  the  peak  in  S/N  changes 

from  5  Hz  at  12°  to  2  Hz  at  23°  (and  continues  at  2  Hz  out  to  33°  as  seen 

from  the  remaining  spectra  in  Appendix  II) ,  the  spectra  in  Appendix  II  also 

show  that  at  NORSAR  10  Hz  energy  is  above  the  noise  all  the  way  out  to  33°. 

The  relative  usefulness  of  the  10  Hz  energy  for  detection  can  only  increase 

for  smaller  events  which  have  higher  corner  frequencies.  Assuming  an  asymptotic 

-2 

spectrum  proportional  to  w  and  a  typical  corner  frequency  near  5  Hz  for  the 
events  studied  (SALMON'S  corner  frequency  was  ~  3  Hz),  then  one  could  regain 
about  12  dB  S/N  at  10  Hz.  This  would  still  not  be  a  better  S/N  than  we  observe 
in  the  range  2  to  5  Hz;  however,  the  noise  at  10  Hz  in  the  spectra  presented  may 
well  be  system  or  quantization  noise,  so  that  the  true  S/N  at  10  Hz  may  be  greater 
than  is  apparent.  If  so,  then  it  is  possible  that  for  a  better  designed  system 
the  S/N  might  be  higher  at  10  Hz  than  at  2  to  5  Hz  for  small  events. 
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SUMMARY 


Amplitude-distance  curves  are  different  in  the  EUS  and  WUS:  P  decays 

-2  5  -3.0  maX 

as  r  '  and  r  '  in  EUS  and  WUS,  while  for  the  maximum  after  3.6  km/sec  on 

-2  -3 

the  vertical  component  (termed  L  )  the  decay  rates  are  r  and  r  .  The  EUS 
results  are  in  general  agreement  with  the  literature  and  with  the  data  presented 
by  Nersesov  and  Rautian  (1964)  for  events  on  the  northern  margin  of  tectonic 
regions  in  the  Soviet  Union  suggesting  that  discrimination  results  in  the  EUS 
are  relevant  to  NSS  stations  within  the  Soviet  Union. 


Using  these  distance  amplitude  relations,  network  mean  amplitudes  at 

1000  km  were  computed  for  and  for  earthquakes  and  explosions  in  the 

EUS  and  WUS  and  a  separation  of  0.6  magnitude  units  was  observed  thus  forming 

a  regional  discriminant.  The  L  is  larger  for  earthquakes  than  for  explosions. 

8 

This  conclusion  is  somewhat  uncertain  in  the  WUS  because  of  the  large  scatter 
in  the  explosion  population.  However,  a  reasonable  explanation  for  this  large 
scatter  is  that  the  small  events  at  NTS  are  at  such  shallow  depths  in  dry 
alluvium  that  the  medium  is  weak,  resulting  in  a  low  corner  frequency.  This 


decreases  the  ratio  P  /L  since  P  contains  comparatively  more  high  frequency 

8  8  8 

than  L^.  The  scatter  probably  would  not  be  a  problem  in  a  true  test  ban  situ¬ 
ation  since  shots  will  be  well  buried  to  avoid  surface  collapse. 


The  WUS  earthquake  P^y  versus  is  displaced  about  0.2  magnitude  units 


The  GNOME  explosion  and  the  Hebgen  Lake  earthquake,  which  are  on  the 
border  between  the  EUS  and  WUS  defined  by  Der,  Massfi  and  Gurski  (1975),  show 
differences  in  amplitude  distance  relations  for  the  same  event  in  different 
provinces.  The  greater  WUS  attenuation  results  in  amplitudes  at  10°  of  about 
0.4  magnitude  units  below  that  in  EUS. 


Analysis  of  the  SALMON  and  18  February  1964  Alabama  earthquake  shows  that 
there  is  no  earthquake/explosion  discrimination  capability  using  maximum  trans¬ 
verse  to  maximum  radial  amplitude  ratios.  We  also  find  that  the  source 


spectra  of  the  two  events  are  identical  between  1  and  10  Hz,  that  the  L^ 
spectrum  is  different  from  the  P  spectrum  and  is  therefore  not  the  source 


spectrum,  and  that  the  L^  spectrum  is  contaminated  by  scattered  coda  from 
earlier  phases  so  that  high-frequencies  observed  in  the  L^  phase  may  not  be 
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predictable  by  any  deterministic  theory  of  L^.  Identical  conclusions  with 
respect  to  the  P  and  spectra  were  obtained  by  analysis  of  spectra  of 
Soviet  explosions  as  observed  at  NORSAR. 

In  the  EUS,  10  Hz  energy  is  observed  out  to  10°,  and  for  Soviet  shots 
observed  at  NORSAR,  out  to  33°.  The  peak  in  the  S/N  varies  from  5  Hz  at 
11°  to  2  Hz  at  33°;  however,  for  small  events  which  have  not  yet  been  studied, 
the  peak  may  well  be  at  higher  frequencies  because  of  higher  corner  frequen¬ 
cies  and  because  much  of  the  noise  on  LRSM  and  NORSAR  systems  at  10  Hz  may 
be  system  or  quantization  noise  which  can  be  reduced  by  more  carefully 
designed  systems. 

Detection  thresholds  for  a  single  element  at  the  C3  subarray  of  NORSAR 
were  determined  to  be  m^  0.9,  2.4,  3.3,  and  3.9  at  5°,  10°,  15°  and  20°, 
respectively.  It  should  be  noted  that  the  lower  frequency  noise  levels  are 
high  at  NORSAR,  which  is  near  the  sea,  and  thus  that  stations  on  continents 
might  have  lower  thresholds  at  large  distances  where  the  lower  frequencies 
are  the  most  useful  for  detection. 
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